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ABSTRACT 


Despite  decades  of  energetic  international  control  efforts,  nuclear  weapons  tech¬ 
nology  continues  to  spread  worldwide.  To  understand  how  these  complex  weapons 
programs  can  be  developed,  we  assume  the  role  of  a  nation  seeking  to  build  a  first  fission 
weapon,  and  the  ability  to  continue  to  build  more.  We  introduce  a  large-scale  project 
management  model  that  includes  alternate  development  paths  to  achieve  certain  key 
technical  milestones.  We  show  how  such  a  project  can  be  optimally  accelerated  by  expe¬ 
diting  critical  tasks.  Next,  we  present  a  new  analysis  tool  to  detect  vulnerabilities  in  such 
a  development  program:  we  seek  optimal  actions  to  impede,  set  back  and/or  otherwise 
frustrate  completion  of  a  first  weapon,  even  if  the  proliferator  knows  what  we  are  doing 
to  delay  things.  This  two-sided  project  evaluation  tool  is  implemented  with  a  combina¬ 
tion  of  off-the-shelf  project  management  software,  optimization  software  and  custom 
code.  An  illustrative  case  study  of  a  first  fission  weapon  program  shows  how  this  new 
analysis  tool  can  be  used.  Our  methods  also  apply  to  chemical,  biological  and/or  radio¬ 
logical  dispersion  weapons,  as  well  as  to  more  conventional  strategic  industrial  and 
commercial  activities. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  the  computer  programs  developed  in  this  research 
may  not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made, 
within  the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic 
errors,  they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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EXECUTIVE  SUMMARY 


We  introduce  a  model  for  planning,  even  expediting,  development  of  a  nation’s 
first  operational  nuclear  weapon.  Our  goal  is  to  understand  how  such  a  complex  project 
can  be  conducted  by  a  belligerent  state  so  that  we  might  find  vulnerabilities  to  exploit: 
we  seek  actions  to  delay,  impede,  interdict  and/or  otherwise  frustrate  completion  of  such 
a  weapon  of  mass  destruction. 

Preventing  the  proliferation  of  such  devices  is  an  international  priority.  The  1968 
Treaty  on  the  Non-proliferation  of  Nuclear  Weapons  (NPT)  has  been  ratified  by  188  na¬ 
tions,  more  than  any  other  international  arms  control  agreement.  Global  proliferation 
concerns  have  been  renewed  by  recent  world  events  such  as  North  Korea’s  withdrawal 
from  NPT  and  revelations  of  Libya’s  clandestine  work  on  nuclear  weapons. 

Nuclear  weapon  development  attracts  intense  global  scrutiny.  The  International 
Atomic  Energy  Agency  (IAEA),  chartered  by  the  United  Nations,  is  responsible  for  veri¬ 
fying  that  nuclear  materials  are  used  only  for  peaceful  purposes  by  enforcing  safeguards 
imposed  by  the  treaty.  Any  nation  that  has  ratified  NPT  and  is  now  attempting  to  de¬ 
velop  its  first  nuclear  weapon  would  therefore  almost  certainly  run  a  covert  program. 
Outside  discovery  of  such  a  program  would  depend  on  intelligence  gathering  or  a  public 
announcement  by  the  proliferators  intended  to  intimidate  other  nations  into  favorable  po¬ 
litical  or  economic  concessions.  When  IAEA  inspections  are  denied  or  a  nation  like 
North  Korea  withdraws  from  the  treaty,  the  international  community  may  take  other  steps 
to  counter  the  proliferation  threat. 

We  introduce  a  project  management  tool  to  identify  the  tasks  in  a  weapons  devel¬ 
opment  program  that  are  most  critical  to  its  timely  completion.  We  highlight  the  effect  of 
delaying  an  individual  task,  or  set  of  tasks,  on  completion  of  a  weapon.  These  delays 
might  be  inflicted  by  deliberate  interdiction  of  the  weapon  project;  in  fact,  we  might 
choose  to  set  back  certain  milestones  already  achieved,  as  Israel  did  in  1981  when  they 
successfully  bombed  Iraq’s  Osirak  reactor  facility.  In  addition  to  military  action,  delay- 
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ing  actions  might  also  include  sabotage  of  facilities  or  supplies,  embargo  or  bloekade  of 
material  shipments,  or  motivating  key  personnel  to  leave  the  program. 

Mathematics-based  methods  of  seheduling  and  managing  eomplex  projeets  have 
been  widely  employed  in  government  and  industry  sinee  the  late  1950’s.  The  original 
Program  Evaluation  Review  Technique/Critieal  Path  Method  (PERT/CPM)  models  have 
been  extended  in  a  variety  of  ways  over  the  years  to  incorporate  situations  that  arise  in 
different  types  of  projeets.  Here,  we  must  combine  several  of  these  individual  embel¬ 
lishments. 

The  weapon  developer,  referred  to  heneeforth  as  the  “proliferator,”  wants  to 
minimize  the  duration  of  the  projeet  to  develop  a  nuelear  weapon  by  redueing  the  dura¬ 
tion  of  as  many  tasks  as  he  ean  afford.  We  present  a  novel  and  flexible  model  that 
eaptures  many  of  the  issues  involved  in  expediting  sueh  a  sophistieated  researeh  and  de¬ 
velopment  projeet.  We  inelude  the  ability  to  ehoose  among  alternatives  that  are 
sometimes  available  for  aehieving  eertain  milestones,  seheduling  restrietions  based  on  the 
availability  of  searee  resourees,  and  “erashing”  the  projeet  (i.e.,  shortening  projeet  dura¬ 
tion  by  applying  additional  resources  to  seleeted  tasks),  among  others.  Elsing  this  tool  as 
a  forecast  of  proliferation  plans,  we  next  seek  an  optimal  set  of  tasks  in  the  project  to  in¬ 
terdict  in  order  to  infliet  the  most  delay  in  its  eompletion. 

The  projeet  interdietor  wants  to  impede  progress  to  make  the  projeet  take  as  long 
as  possible,  subjeet  to  eonstraints  on  what  he  can  afford  to  do.  He  may  exert  interdietion 
effort  against  a  partieular  task  in  a  variety  of  ways,  eonstrained  by  a  monetary  budget  or 
the  politieal,  environmental  or  eeonomie  implieations  of  his  aetions. 

We  develop  a  two-sided  optimization:  the  proliferator  and  the  interdietor  can 
each  see  what  the  opponent  is  doing,  and  ean  eaeh  behave  optimally  with  resourees  at 
hand  to,  respeetively,  minimize  and  maximize  the  time  to  eompletion  of  a  first  fission 
weapon.  This  optimization  provides  a  powerful  analytieal  tool.  It  ean  be  used  to  identify 
the  most  vulnerable  tasks  in  the  project,  to  suggest  alternate  delay  plans  with  equivalent 
delaying  effects  to  a  policy  maker  and  to  explore  the  sensitivity  of  the  projeet  duration  to 
the  inputs. 
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Our  prototypic  decision  support  system  is  implemented  using  off-the-shelf  soft¬ 
ware.  All  user  interaction  and  data  management  are  done  through  software  already 
familiar  to  millions  of  personal  computer  users.  Powerful  optimization  software,  unfa¬ 
miliar  to  many  potential  users  of  this  product,  is  hidden  behind  a  straightforward  user 
interface  incorporated  within  standard  project  management  software. 

Computational  experience  with  our  case  study  shows  that  good  solutions  can  be 
found  quickly.  “What-if  ’  questions  regarding  particular  delay  plans  can  be  answered  in  a 
fraction  of  a  second  using  a  desktop  personal  computer.  An  analyst  can  take  advantage 
of  this  speed  to  explore  a  wide  range  of  assumptions  and  find  the  tasks  that  are  most  criti¬ 
cal  to  the  expeditious  completion  of  the  proliferator’s  project. 

This  low-cost  tool  can  be  used  today.  It  provides  quantitative  information  that 
can  be  used  to  guide  policy  makers  working  to  stop  the  proliferation  of  nuclear,  chemical, 
biological  and/or  radiological  dispersion  weapons,  or  to  disrupt  more  conventional  strate¬ 
gic  industrial  and  commercial  activities. 

In  our  case  study,  an  uninterdicted  first  weapon  is  completed  in  260  weeks  (5 
years).  Interdiction  of  this  proliferation  program,  if  kept  secret  from  the  proliferator,  de¬ 
lays  completion  to  340  weeks  (6X  years).  Finally,  if  we  assume  that  the  proliferator  can 
see  what  the  interdictor  is  doing,  and  adjust  his  plan  accordingly,  we  still  delay  comple¬ 
tion  to  324.4  weeks  (6X  years).  This  demonstrates  the  value  to  the  interdictor  of 
keeping  his  actions  secret,  but,  more  importantly,  exposes  the  inherent  vulnerability  of 
such  a  large  project  to  a  few  well-placed  interdiction  actions.  It  also  demonstrates  to  the 
proliferator  that  a  covert  project  can  be  completed  nearly  15  months  faster  than  one 
which  is  known  by  the  enemy  and  therefore  subject  to  interdiction. 

We  are  suggesting  how  to  organize  and  manage  intelligence  data  in  a  cohesive 
fashion  that  offers  visualization,  data  base  functionality  and  interoperability  with  Micro¬ 
soft  Office.  The  optimization  is  new,  but  hidden  cleanly  from  view.  An  analyst,  or  even 
a  policy  maker,  can  pose  questions,  get  answers  and  gain  deep  insights  without  any 
mathematical  background  at  all. 
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I.  NUCLEAR  PROLIFERATION 


A.  BELLIGERENT  STATES  SEEK  NUCLEAR  WEAPONS 

We  introduce  a  model  for  planning,  even  expediting,  development  of  a  nuclear 
weapon.  Our  goal  is  to  understand  how  such  a  complex  project  can  be  conducted  by  a 
belligerent  state  so  that  we  might  find  vulnerabilities  to  exploit:  we  seek  actions  to  delay, 
impede,  interdict  and/or  otherwise  frustrate  completion  of  such  a  weapon  of  mass  de¬ 
struction. 

Nearly  sixty  years  after  the  first  nuclear  weapon  was  detonated,  preventing  the 
proliferation  of  such  devices  is  an  international  priority.  One  hundred  eighty-eight  na¬ 
tions  have  ratified  the  1968  Treaty  on  the  Non-proliferation  of  Nuclear  Weapons  (NPT), 
more  than  any  other  international  arms  control  agreement  [United  Nations,  2004].  North 
Korea’s  recent  withdrawal  from  NPT  and  Libya’s  revelations  of  work  on  nuclear  weap¬ 
ons  have  renewed  global  proliferation  concerns. 

Production  of  nuclear  weapons  attracts  intense  global  scrutiny;  because  of  NPT 
safeguards  currently  in  place,  any  nation  attempting  to  develop  its  first  nuclear  weapon 
would  almost  certainly  run  a  covert  program.  Outside  discovery  of  such  a  program 
would  depend  on  intelligence  gathering  or  a  public  announcement  by  the  proliferators 
intended  to  intimidate  other  nations  into  favorable  political  or  economic  concessions. 
The  International  Atomic  Energy  Agency  (IAEA),  chartered  by  the  United  Nations,  is 
responsible  for  verifying  that  nuclear  materials  are  used  only  for  peaceful  purposes 
[United  Nations,  1968,  Article  III].  When  IAEA  inspections  are  denied  or  a  nation  like 
North  Korea  withdraws  from  the  treaty,  the  international  community  may  take  other 
steps:  the  successful  attack  on  Iraq’s  Osirak  nuclear  facility  by  the  Israeli  Air  Eorce  in 
1981  [e.g.,  Eederation  of  American  Scientists,  2000]  is  an  example.  Alarmed  by  intelli¬ 
gence  that  Iraq  could  assemble  a  functional  nuclear  weapon  within  two  years,  the  Israelis 
chose  an  immediate  military  option  rather  than  a  diplomatic  solution. 

B,  A  PROPOSED  COURSE  OF  ACTION 

We  introduce  a  project  management  tool  to  identify  the  tasks  in  a  weapons  devel¬ 
opment  program  that  are  most  critical  to  its  timely  completion.  We  highlight  the  effect  of 
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delaying  an  individual  task,  or  set  of  tasks,  on  eompletion  of  a  weapon.  These  delays 
might  be  inflieted  by  deliberate  interdietion  of  the  weapon  projeet:  in  faet,  we  might 
choose  to  set  back  certain  milestones  already  achieved.  Delaying  aetions  might  inelude 
sabotage  of  faeilities  or  supplies,  embargo  or  bloekade  of  material  shipments,  or  even 
military  aets.  Although  our  example  is  nuelear  weapon  development,  our  methods  apply 
just  as  well  to  ehemieal,  biologieal  and/or  radiologieal  dispersion  weapons,  as  well  as  to 
more  eonventional  strategie  industrial  and  eommereial  aetivities. 

Sinee  the  late  1950’s,  mathematies-based  methods  of  seheduling  and  managing 
eomplex  projeets  have  been  widely  employed  in  government  and  industry.  The  original 
Program  Evaluation  Review  Teehnique/Critieal  Path  Method  (PERT/CPM)  models  [Mal- 
eolm,  et  al,  1959]  have  been  extended  over  the  years  to  ineorporate  a  variety  of  situations 
that  arise  in  different  types  of  projeets,  sueh  as  alternate  ehoiees  that  are  sometimes  avail¬ 
able  for  aehieving  eertain  milestones,  seheduling  restrietions  based  on  the  availability  of 
searee  resourees,  and  “erashing”  the  projeet  (i.e.,  shortening  projeet  duration  by  applying 
additional  effort  to  seleeted  tasks),  among  other  additions.  Due  to  the  eomplexity  of  nu- 
elear  weapons  programs,  we  need  all  these  embellishments. 

Projeet  management  models  are  universally  represented  as  networks  [e.g.,  Moder, 
Phillips  and  Davis,  1983,  Chapter  1].  In  sueh  a  representation,  the  overall  duration  of  the 
projeet  is  the  additive  length  of  the  longest  path  through  the  network,  also  known  as  the 
critical  path.  The  weapon  developer,  referred  to  heneeforth  as  the  “proliferator,”  wants 
to  minimize  the  duration  of  the  projeet  to  develop  a  nuelear  weapon  by  redueing  the  dura¬ 
tion  of  as  many  “eritieal”  tasks  as  he  ean  afford.  While  building  the  projeet,  he  may  be 
able  to  ehoose  among  several  alternate  means  available  to  reaeh  some  milestone. 

The  projeet  interdietor  wants  to  impede  progress  to  make  the  projeet  take  as  long 
as  possible,  subjeet  to  eonstraints  on  what  the  interdietor  ean  afford  to  do.  He  may  exert 
interdietion  effort  against  a  partieular  task  in  a  variety  of  ways,  eonstrained  by  a  mone¬ 
tary  budget  or  the  politieal,  environmental  or  eeonomie  implieations  of  his  aetions. 

We  develop  a  two-sided  optimization;  the  proliferator  and  the  interdietor  ean 
eaeh  see  what  their  opponent  is  doing,  and  ean  eaeh  behave  optimally  with  resourees  at 
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hand  to,  respectively,  minimize  and  maximize  the  completion  of  a  first  fission  weapon, 
and  the  ability  to  continue  producing  more. 
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II.  MODELING  PROLIFERATION  AND  ITS  INTERDICTION 


A.  PROJECT  MANAGEMENT  IN  NUCLEAR  PROLIFERATION 

Embarking  on  a  program  to  develop  a  nation’s  first  nuelear  weapon  requires  a 
signifieant  investment  of  materials,  manpower  and  technology.  Even  for  a  modern  indus¬ 
trialized  nation,  this  investment  will  likely  represent  a  significant  fraction  of  that  nation’s 
available  resources.  Such  a  program  is  a  massive  undertaking  and  we  assume  that  it 
would  be  subject  to  intense  high-level  management  scrutiny  and  careful,  centralized  co¬ 
ordination. 

Moder,  Phillips  and  Davis  [1983,  p.  3]  define  a  project  as  “a  set  of  tasks  or  activi¬ 
ties  related  to  the  achievement  of  some  planned  objective,  normally  where  the  objective 
is  unique  or  non-repetitive.”  Our  proliferator’s  program  to  develop  his  first  nuclear 
weapon  fits  this  definition  well.  We  can  reasonably  expect  him  to  employ  standard  pro¬ 
ject  management  tools  to  plan  and  organize  the  project  in  order  to  efficiently  schedule  its 
tasks.  We  assume  that  once  the  proliferator  commits  substantial  national  resources  to¬ 
wards  the  development  of  a  nuclear  weapon,  he  wants  to  complete  his  project  and  obtain 
the  first  operational  weapon  as  quickly  as  possible. 

Our  proliferator’s  problem  can  be  represented  most  simply  as  a  classic  project 
network  enhanced  to  allow  crashing  [e.g.,  Moder,  Phillips  and  Davis,  1983,  Chapter  8]. 
However,  we  recognize  that  this  model  is  inadequate  to  represent  many  situations  that 
can  occur  in  practical  real-world  projects.  To  this  classic  project  network,  we  add  the 
following  embellishments: 

•  Completion  of  any  task  in  a  “normal”  amount  of  time  consumes  a  fixed 
amount  of  one  or  more  non-renewable  resources. 

•  The  duration  of  an  individual  task  may  be  shortened  (“crashed”)  by  allocation 
of  additional  quantities  of  the  required  resources.  We  assume  a  linear  rela¬ 
tionship  between  the  amounts  of  additional  resources  provided  and  the 
duration  of  the  task — more  resources  accelerates  progress.  Each  task  requires 
some  minimum  amount  of  time,  no  matter  how  many  additional  resources  are 
allocated  to  it. 

•  Crashing  is  limited  by  the  availability  of  each  resource  and  by  an  overall 
monetary  budget. 
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•  Certain  milestone  tasks  may  be  aehieved  via  alternate  eourses  of  aetion. 
When  one  alternate  path  is  ehosen,  the  tasks  in  the  other  alternative(s)  need 
not  be  eompleted.  A  decision  node  is  a  task  in  the  projeet  network  where  al¬ 
ternate  eourses  of  aetion  diverge. 

•  The  project  being  delayed  may  already  be  partially  completed.  If  a  task  that 
was  previously  in  progress  or  completed  is  interdicted,  it  and  all  of  its  succes¬ 
sors  must  be  restarted  and  completed  again  in  order  to  finish  the  project. 

•  Standard  finish-to-start  precedence  relationships  between  pairs  of  tasks  are 
generalized  to  include  start-to-start,  finish-to-finish  or  start-to-finish.  There  is 
a  finite  time  lag  (or  lead)  associated  with  each  precedence  constraint. 

•  The  project  may  include  a  hierarchical  grouping  of  the  tasks.  A  summary  task 
represents  the  tasks  grouped  below  it  in  a  hierarchy  (its  “children”)  but  does 
not  consume  resources  or  otherwise  contribute  to  the  completion  of  the  pro¬ 
ject.  Summary  tasks  can  be  given  their  own  precedence  relationship(s)  with 
other  tasks  in  the  project.  Accommodating  them  requires  introducing  dummy 
precedence  arcs  into  the  network  to  ensure  a  summary  task  starts  no  later  than 
any  of  its  children  and  to  preclude  finishing  until  all  of  its  children  finish.  The 
summary  task  must  be  completed  if  any  of  its  children  are  completed  and  the 
children  must  be  completed  if  the  summary  task  is  forced  to  start  by  other 
precedence  constraints. 

•  Interdicting  a  task  delays  its  completion  by  a  constant  amount. 

•  Possible  delay  plans  are  limited  by  a  maximum  total  number  of  interdictions 
and  by  an  interdiction  “budget.”  This  budget  may  represent  total  financial 
cost  or  level  of  geopolitical  impact  for  each  interdiction  action. 

The  proliferator’s  objective  is  to  minimize  the  length  of  the  critical  path  in  this 
network.  We  first  determine  what  the  proliferator  would  do  without  any  threat  of  inter¬ 
diction. 

B,  EXPEDITING  A  PROJECT  WITH  NO  INTERDICTION 

We  are  now  going  to  formalize  these  concepts  with  a  generalization  of  a  well- 
known  project  management  tool. 

1.  Sets  and  Index  Use 

We  use  an  activity  on  node  (AON)  formulation  of  the  project  network,  employing 
the  following  sets  of  tasks,  task  relationships  and  resources.  In  the  AON  formulation,  the 
nodes  in  the  network  represent  the  individual  tasks  in  the  project  and  the  arcs  represent 
the  pairwise  partial  orders  between  respective  predecessor  and  successor  tasks. 
i,j  e  N  tasks  (nodes) 

cr,  r  e  A  start  and  finish  milestone  tasks,  respectively 
DN  Cl  N  decision  nodes 
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STc^N 

AN^N 

(U)  e  A 

FS^A 

FF^A 

SF^A 

SS^A 

CHILD  ^  A 
r^R 


hierarchical  summary  tasks 

all  tasks  except  the  decision  nodes  and  summary  tasks; 
AN  =  N-DN-ST 

precedence  relationships:  task  i  precedes  task  j  (arcs) 

(ij')  e  FS  implies  task  i  must  finish  before  task  j  can  start 

(ij')  e  FF  implies  task  i  must  finish  before  task  j  can  finish 

(ij')  e  SF  implies  task  i  must  start  before  task  j  can  finish 

(ij')  e  SS  implies  task  i  must  start  before  task  j  can  start 

FS,  FF,  SF  and  SS  are  mutually  exclusive  and  exhaustive  subsets  of  A 

(ij')  e  CHILD  if  task  j  is  subordinate  to  summary  task  i 

non-renewable  resources 


2.  Data  [units] 


t 

Pi 

lagij 
d.,  J, 


b,  b, 


rr 

LI—u 


rer. 


Qr 

budget 

delayi 

m 


3. 


elapsed  time  since  the  beginning  of  the  project  [days] 
fraction  of  work  on  task  i  ^  N  that  is  complete  at  time  t  [  0 
time  lag  for  a  precedence  constraint  e  A  [days] 


duration  of  task  i  ^  N  with  no  crashing  effort  and  maximal  crashing  ef¬ 
fort,  respectively  [days] 

minimum  and  maximum  number,  respectively,  of  successor  tasks  allowed 
to  be  chosen  at  decision  node  i  e  DN  [non-negative  integer] 
minimum  resource  requirement  for  resource  r  e  i?  to  complete  task  i  ^  N 
with  no  crashing  effort  [resource  r-units] 

resource  expediting  rate  of  resource  r  e  i?  for  task  i  ^  N  [resource  r- 
units/day] 

unit  cost  of  resource  r  e  i?  [$/resource  r-unit] 
total  availability  of  resource  r  e  i?  [resource  r-units] 
crashing  cost  budget  [$] 

amount  task  i  e  N  -  DN  -  ST  -  <j-  t  is  delayed  when  interdicted  [days] 
a  “large”  number  used  to  relax  the  precedence  constraints  for  tasks  that  are 

not  actually  completed;  m  =  maxCd^  +  delay max  )  [days] 


Variables  [units] 


S.^  earliest  start  time  of  task  i  e  N  [days] 

D.  the  decision  plan;  Z).  =  1  if  task  /  e  V  is  completed,  0  otherwise  [binary] 

4  the  crashing  plan;  5^  is  the  amount  task  i  e  N -DN -<j-  r  is  accelerated 

by  crashing  [days] 
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4,  Expediting  a  Project 

min  5,- (aO) 
s.t.  Sj  +  m(l-D^)>S.  +  (d.-S.'J+lagy  ^(ij)eFS  (al) 

S,+^,-Sj)+m(l-D,)>  y(i,j)^FF  (a2) 

5,  +  (d,-S^+lag,j 

Sj  +  (d  .  -  m (1  -  D.)  >  S.  +  lag.j  V (/ j')  e  SF  (a3) 

5^+m(l-Z),)>5,  +  /ag,.  V(/j)e55  (a4) 

[}Ti,.D.  +  rer-^.S. )  <  budget  (a5) 

r  L  ieN-<7—T 

'^{iZirDi  +  rer^^d^<  Vr  e  R  (a6) 

i^N-a  —  T 

bp.  <  Y,Dj  <  bp,  Vz  e  DN  (a7) 

j  I  (i,j)£A-(CHlLDnFF) 

D.  >  D,  V(z  J-)  ^  A-  CHILD  \i  ^  {N  -  DN)  (a8) 

Dj  >  D,  V(z  J-)  e  CHILD  \  (a9) 

'(3(zz,z)  e  A  I  ((zz,z)  ^  CHILD &{i,n)  ^  CHILD)) 

^&(not  S(nj)  e  A  \  n  e  DN)  ^ 

0<Spd.-di  \/i  ^  N-DN-<j- T  (alO) 

Z),  e  {0,1}  yteN  (all) 

5^=0,Z)^=1  (al2) 

D,=\  yieN\p,>0  (al3) 

Spt  yt  e  N-ST\p.  =  0  (al4) 

Si  =  t- Pi(d^  -Si)  Vz  e  1 0  <  />,  <  1  (al5) 

Si<t-(di-d^  'yi^N\pi=\  (al6) 


The  proliferator  seeks  to  minimize  the  overall  length  of  the  projeet,  expressed  in 
(aO)  as  the  differenee  between  the  finish  time  and  the  start  time.  All  of  the  choiees  made 
among  alternatives  at  the  deeision  nodes  are  eaptured  in  the  resulting  decision  plan.  The 
crashing  plan  speeifies  the  level  of  additional  effort  applied  to  eaeh  task  to  speed  it  up. 
The  deeision  plan  and  the  erashing  plan  together  eneompass  all  of  the  deeisions  that  the 
proliferator  must  make  to  expedite  the  projeet.  They  are  eolleetively  referred  to  as  the 
expediting  plan. 

Constraints  (al)  to  (a4)  enforee  pairwise  task  partial  orders.  Eaeh  of  these  eon- 
straints  features  the  length  of  a  task,  -  Si),  its  normal  duration  minus  the  amount  of 

erashing  effort.  The  term  zn(l-Z).)  [Crowston  and  Thompson,  1967]  relaxes  its  eon- 
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straint  if  task  i  is  not  one  of  the  tasks  aetually  eompleted  in  the  final  deeision  plan  (i.e.,  a 
deeision  takes  an  alternate  eourse  of  aetion). 

Constraints  (a5)  and  (a6)  respeetively  enforee  the  monetary  budget  and  resouree 
availabilities.  Constraints  (a7)  limit  the  number  of  alternatives  allowed  at  eaeh  deeision 
point.  Constraints  (a8)  ensure  that  a  task  is  foreed  to  be  eompleted  if  any  of  its  predeees- 
sors  are,  exeept  for  tasks  that  are  ehildren  of  a  summary  task.  Constraints  (a9)  foree 
ehildren  of  a  summary  task  to  be  eompleted  only  if  the  summary  task  has  predeeessors 
and  the  ehild  is  not  direetly  affeeted  by  a  deeision.  (a  10)  show  that  eaeh  task  has  a 
maximum  erash  effort  (that  gives  it  a  shortest  duration).  Task  eompletion  deeisions  are 
required  to  be  binary  in  (all).  Constraint  (al2)  eompletes  the  starting  task  no  matter 
what  deeisions  are  made  later  on  and  starts  the  projeet  at  time  zero. 

Constraints  (al3)  foree  the  deeision  plan  to  inelude  every  task  that  has  already 
eommeneed,  re  fleeting  any  deeisions  already  made.  Constraints  (a  14)  require  tasks  with 
no  reeorded  progress  to  start  no  sooner  than  the  eurrent  time.  (al5)  set  the  start  time  of 
tasks  that  have  started  but  not  finished  while  (al6)  ensure  eompleted  tasks  are  assigned 
start  times  early  enough  to  have  finished  by  the  eurrent  time. 

5,  Model  Derivation 

Eaeh  task  in  the  projeet  has  a  normal  duration,  the  amount  of  time  required  to 
eomplete  without  any  speeial  expediting  effort.  Completing  a  task  requires  the  expendi¬ 
ture  of  some  fixed  amount  of  one  or  more  non-renewable  resourees  (raw  materials, 
energy,  labor,  maehinery,  ete.).  A  fixed  amount  of  eaeh  resouree  is  alloeated  to  the  entire 
projeet.  The  duration  of  a  task  ean  be  shortened  down  to  some  non-negative  minimum 
level  by  applying  additional  resourees  to  “erash”  the  task.  The  proliferator  may  expend 
his  resourees  at  any  time  and  any  rate  he  ehooses,  but  the  eumulative  usage  of  eaeh  re- 
souree  may  never  exeeed  the  total  quantity  available.  If  the  projeet  is  interdieted  and  a 
task  is  delayed,  a  fixed  amount  of  time  is  added  to  the  duration  of  that  task.  The  relation¬ 
ship  between  these  times  is  shown  in  Figure  1. 
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Expediting  effort 


Task  start 


1 

Minimum 
task  duration 

Maxi 

task  dt 

mum 

iration 

Task  finish 


Interdiction  delay  Expediting  effort 


1  F 

* 

Minimum 

Task  start  task  duration  Task  finish 


Figure  1.  Task  duration  with  expediting  effort,  both  without  and  including  inter¬ 
diction  delay. 

The  duration  of  an  uninterdicted  task  {left)  may  be  expedited  from  its  longest 
duration  to  its  shortest.  Interdicting  the  task  adds  a  fixed  length  of  time  to  its 
duration  {right)  while  still  allowing  expediting  effort  to  reduce  its  overall  du¬ 
ration. 


Early  analytical  project  management  work  [e.g.,  Malcolm,  et  al,  1959]  assumed 
that  if  a  task  depended  on  another  in  any  way,  the  later  task  could  not  begin  until  imme¬ 
diately  after  all  of  its  prerequisites  had  been  completed.  This  was  later  generalized  in  a 
model  called  precedence  diagramming  by  J.  David  Craig  [International  Business  Ma¬ 
chines,  1968].  Precedence  diagramming  permits  every  combination  of  pairwise  partial 
orders  between  predecessor  and  successor  task  start  and  finish.  Given  this  generality,  it  is 
easy  to  accommodate  a  lag-  or  lead-time  between  pairs  of  tasks.  These  relationships  are 
shown  in  Figure  2. 


S  1  F 
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S  2  F 

(a) 

S  1  F 
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(c) 
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Figure  2,  Partial  orders  between  predecessor  and  successor  tasks. 

Possible  combinations  are  finish-to-start  {a),  start-to-start  {b),  finish-to-finish 
(c)  and  start-to -finish  {d).  All  lags  are  shown  as  positive  values  (Task  2  “lags” 
Task  1),  although  zero  and  negative  lags  (“leads”)  are  also  allowed. 


Crowston  and  Thompson  [1967]  introduce  decision  CPM  that  allows  alternate 
courses  of  action  to  achieve  milestones.  Their  formulation  divides  the  project  tasks  into 
sets  based  on  where  the  decision  points  are  in  the  project  network.  Each  set  of  tasks  is 
assigned  a  binary  variable  that  indicates  whether  that  set  of  tasks  is  actually  completed  in 
the  project.  For  our  purposes,  it  is  more  practical  to  assign  each  task  in  the  project  its 
own  binary  decision  variable  and  add  additional  constraints  to  maintain  the  decision  set 
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Figure  3.  A  decision  critical  path  method  project  network  depicting  one  decision. 

Similar  to  the  notation  introduced  by  Crowston  and  Thompson  [1967],  deci¬ 
sion  nodes  are  indicated  by  triangle.  The  nodes  with  the  heavy  border  must  be 
completed  in  order  to  finish  the  project.  If  Task  7  is  chosen  {left),  all  tasks  are 
completed  while  if  Task  8  is  chosen  {right).  Task  7  is  never  completed.  Dot¬ 
ted  lines  indicate  discretionary  precedence  constraints  that  are  relaxed. 

relationships  between  them.  The  effect  of  a  decision  on  the  completion  of  a  project  is 

shown  in  Figure  3. 

C.  DELAYING  A  PROJECT  WITH  A  FIXED  DECISION  PLAN 

The  interdictor  chooses  a  set  of  project  tasks  to  delay  in  order  to  maximize  the 
length  of  the  critical  path.  This  model  assumes  that  the  proliferator  is  unaware  of  the  in¬ 
terdiction  effort  and  does  nothing  to  compensate  even  after  the  interdictor  inflicts  a  delay. 
This  requires  as  input  the  result  of  expediting  the  project  with  no  interdiction.  Slack  is 
calculated  for  each  task  from  an  earliest  start  time  and  a  latest  start  time  as  in  conven¬ 
tional  CPM.  Expediting  the  project  provides  the  earliest  start  times  directly  and  latest 
start  times  are  easily  calculated  using  the  length  of  each  expedited  task.  Tasks  on  the 
critical  path  have  zero  slack  while  a  positive  slack  implies  the  task  is  not  on  the  critical 
path. 
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1. 


Data  [units] 

dCosti  cost  to  delay  task  i  e  AN -a-  t  [$] 

dBudget  delaying  eost  budget  [$] 

maxint  maximum  number  of  interdietions  allowed  [eount] 

slacks  the  differenee  between  the  latest  start  time  and  earliest  start  time  for  task 

i  e  N  [days]. 

d.  the  fixed  deeision  plan;  d.  =  1  if  task  i  ^  N  is  aetually  eompleted,  0  oth¬ 

erwise  [binary] 

2.  Variables  [units] 

X.  the  delay  plan;  X,  =l  if  task  i  e  AN-a-r  is  interdieted,  0  otherwise 


[binary] 

3,  Delaying  a  Project 

max  '^(delay^- slack  (bO) 

i\di  =  l 

<  maxint  (bl) 

iE.AN-<j-T 

'y' ,dCost  -X ;  <  dBudget  (b2) 

iE.AN-a—T 

X.  e  {0,1}  V/  e  AN -a-  t  (b3) 


The  objeetive  funetion  (bO)  greedily  favors  delay  plans  that  he  on  the  eritieal  path 
{d.  =  \  and  slack-  =  0).  The  delay  added  to  the  projeet  is  either  the  total  delay  from  that 
task  if  it  is  on  the  expedited  eritieal  path  or  an  amount  redueed  by  the  slaek  at  that  task. 
Constraint  (bl)  limits  the  maximum  number  of  interdietions  allowed  and  (b2)  enforees 
the  delaying  budget.  Constraints  (b3)  require  the  interdietion  deeisions  to  be  binary. 

Even  though  the  interdietor  may  aehieve  eomplete  surprise  with  his  delaying  ef¬ 
fort,  as  Israel  did  with  its  attaek  on  Iraq’s  Osirak  reaetor  faeility,  the  proliferator  will 
almost  eertainly  become  aware  of  the  interdietion  after  the  faet.  We  assume  an  intelligent 
adversary  will  reevaluate  his  options  at  this  point  and  take  immediate  aetion  to  mitigate 
the  effeets  of  the  interdietion. 

D,  A  TWO-SIDED  MODEL:  DELAYING  A  PROJECT  WHEN  THE 
PROLIFERATOR  AND  THE  INTERDICTOR  CAN  SEE  EACH  OTHER’S 
ACTIONS 

1.  Finding  the  Best  Delay  Plan 

We  now  assume  eomplete  transparency  between  proliferator  and  interdietor — 

eaeh  side  is  fully  aware  of  the  other’s  eapabilities,  limitations  and  intentions.  We  seek  a 

delay  plan  that  maximizes  the  length  of  the  projeet  even  if  the  proliferator  antieipates  our 
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1)  Start  with  proliferator’s  optimal  plan  without  interdiction 

2)  Find  a  unique  interdiction  plan 

3)  Determine  the  proliferator’s  response  to  this  plan 

4)  Repeat  steps  2)  through  4)  as  long  as  desired  or  until  all  possible  delay 
plans  have  been  evaluated 

5)  Report  best  solution  discovered 

Figure  4,  Enumeration  procedure. 

actions  and  responds  optimally.  We  implement  decomposition  methods  to  accomplish 
this — ^we  enumerate  delay  plans  and  then  expedite  the  interdicted  project.  We  refer  to  the 
enumeration  of  delay  plans  as  the  master  problem  and  to  expediting  an  interdicted  project 
as  the  subproblem. 

Figure  4  outlines  our  enumeration  procedure.  We  start  by  determining  the  opti¬ 
mal  starting  times  for  each  task  in  the  project  when  there  is  no  interdiction  effort 
whatsoever.  The  alternatives  chosen,  task  start  times,  total  slack  for  each  task,  resource 
usage  levels,  the  crashing  plan  employed  and  the  length  of  the  project  are  recorded,  estab¬ 
lishing  an  incumbent  baseline  to  use  later.  Information  about  the  task  slack  identifies  the 
critical  path  and  enables  the  master  problem  to  quickly  find  candidate  delay  plans  that 
will  have  a  large  impact  on  the  project  length.  Infeasibility  in  the  subproblem  at  this 
stage  indicates  a  problem  with  the  data — insufficient  budget  or  resources,  conflict  be¬ 
tween  the  elapsed  time  and  progress  of  certain  tasks,  or  et  cetera.  If  the  subproblem  is 
feasible,  we  proceed  by  initializing  a  counter  that  tracks  the  number  of  times  we  solve  the 
master  problem. 

Our  master  problem  finds  a  delay  plan  with  the  maximum  effect  on  the  project 
length  that  is  different  in  at  least  one  detail  from  all  the  delay  plans  that  have  already 
been  considered.  For  a  task  that  is  on  the  critical  path,  delaying  it  adds  the  full  interdic¬ 
tion  delay  to  the  length  of  the  project.  If  a  task  is  not  on  the  critical  path,  an  interdiction 
has  to  lengthen  it  enough  to  overcome  its  slack  duration  before  adding  anything  to  the 
length  of  the  project.  Both  of  these  statements  assume  that  the  proliferator  does  not  then 
choose  a  different  course  of  action  that  excludes  the  task  in  question.  Infeasibility  in  the 
master  problem  indicates  that  we  have  enumerated  every  feasible  delay  plan. 
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The  subproblem  evaluates  the  latest  eandidate  delay  plan  generated  by  the  master 
problem.  If  this  delay  plan  interdiets  all  of  the  tasks  aetually  eompleted  in  the  expedited 
projeet,  or  if  every  interdieted  task  is  on  the  resulting  eritieal  path,  then  we  have  found  a 
lower  bound  on  the  amount  of  delay  that  the  interdietor  ean  infliet.  We  eompute  the 
bound  by  reeomputing  the  master  problem  objeetive  using  the  new  slaek  values  from  the 
subproblem.  We  require  all  subsequent  eandidate  delay  plans  to  eause  at  least  this  mueh 
delay  to  the  projeet.  If  the  subproblem  is  infeasible,  the  eurrent  delay  plan  is  ideal — the 
proliferator  has  exhausted  all  of  his  resources  and  can  never  finish  the  project.  We  re- 
eord  the  length  of  the  delayed  projeet  (infinity  if  it  was  infeasible)  for  this  delay  plan. 
Expeeting  that  there  may  be  alternate  delay  plans  that  lead  to  an  identieal  or  improved 
expedited  projeet  length,  we  eontinue  to  solve  the  master  problem  to  find  another  eandi¬ 
date  delay  plan. 


Eaeh  eyele  of  solving  the  master  problem  followed  by  the  subproblem  is  termed 
an  iteration.  This  eyele  repeats  until  all  interdietion  options  have  been  explored  and  the 
optimal  delay  plan  ean  be  determined,  the  master  problem  objeetive  falls  below  the  iden¬ 
tified  lower  bound,  or  the  number  of  iterations  reaehes  some  preset  limit.  Beeause  we 
save  information  about  the  eandidate  results  of  eaeh  iteration,  it  is  easy  to  find  and  report 
the  best  eandidates. 


2.  Data  [units] 


budgetUsed 

iter 

maxiter 

totalDelay 

lengthk 

slack.^. 


d;, 


erashing  eost  budget  eonsumed  prior  to  time  t  [$] 

eurrent  iteration  eounter  [eount] 

maximum  number  of  iterations  allowed  [eount] 

total  delay  to  projeet  tasks  eaused  by  a  delay  plan  [days] 

the  overall  expedited  projeet  length  in  iteration  k  e  f),..., maxiter]  [days] 

the  differenee  between  the  latest  start  time  and  earliest  start  time  for  task 

i  e  N  in  iteration  k  e  {D,..., maxiter]  [days]. 

amount  of  resouree  r  e  i?  eonsumed  prior  to  time  t  [resouree  r-units] 

the  deeision  plan  for  iteration  k  e  f),..., maxiter]-,  d.^  =1  if  task  i  ^  N  is 

aetually  eompleted  in  iteration  k,  0  otherwise  [binary] 

earliest  start  time  of  task  /  e  in  iteration  k  e  {f,..., maxiter]  [days] 

=  1  if  task  /  e  cr  is  a  sueeessor  of  any  task  in  the  delay  plan  being 
evaluated,  0  otherwise  [binary] 

the  erashing  plan  for  iteration  k  e  maxiter]-,  5.^.  is  the  amount  of 
time  saved  on  task  i  in  iteration  k  by  expediting  its  eompletion  [days] 
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x.^.  the  delay  plan  for  iteration  A:  e  {0,..., max/ter};  x.j^=\  if  task 

i  e  AN  -a  -  t  is  interdieted  in  iteration  k,  0  otherwise  [binary] 

3,  Variables  [units] 

f/.  f/,  =  1  if  task  /  e  V  -  cr  is  in  the  eurrent  delay  plan  or  is  a  sueeessor  of  any 

task  in  the  delay  plan,  0  otherwise  [binary] 

4,  Enumerating  Delay  Plans 

The  interdietor  ehooses  a  set  of  projeet  tasks  to  delay  to  maximize  the  length  of 
the  longest  path  in  the  projeet  network.  This  model  finds  a  set  of  tasks  to  delay. 

max  y^idelay.  -  slack,  +  '^delay.X-  (eO) 

A^ijter-X  ^i,iter-\  ~  ^ 

SM  Z(l-Jf,)+  Xa-,>1  Vk  e  {0,...,/ter-l}  (el) 

i\x,t=\  i\Xit  =  0 

<  maxint  (e2) 

iE.AN-(j-T 

'y' ,dCost ,X ,  <  dBud^et  (e3) 

iEiAN-a-T 

X.  e  {0,1}  Vz  e  AN -a-  r  (e4) 

Given  the  newest  expediting  plan,  the  objeetive  function  (eO)  evaluates  candidate 
delay  plans.  If  a  task  is  completed  in  the  newest  expedited  plan,  then  the  new  delay 
added  is  either  the  total  delay  from  that  task  if  it  was  on  the  expedited  critical  path  or  an 
amount  reduced  by  the  slack  at  that  task.  Tasks  not  included  in  the  newest  expediting 
plan  are  given  a  maximal  delay  because,  in  some  situations,  a  combination  of  other  inter¬ 
dictions  could  force  the  proliferator  to  choose  a  different  decision  plan  that  includes  tasks 
that  were  previously  excluded. 

There  is  a  diversity  constraint  (el)  for  each  legacy  delay  plan  that  requires  the 
Hamming  distance  [Hamming,  1986]  between  each  legacy  delay  plan  and  any  new  delay 
plan  to  be  at  least  one  (i.e.,  a  new  delay  plan  must  differ  in  at  least  one  detail).  Master 
problems  with  diversity  constraints  are  not  new  [e.g.,  see  the  “covering  decomposition” 
of  Israeli  and  Wood,  2002].  Constraints  (e2)  through  (e4)  are  identical  to  constraints  (bl) 
through  (b3). 

Our  enumeration  algorithm  will  eventually  evaluate  every  feasible  delay  plan  in 
the  search  for  an  optimal  solution.  We  stop  the  enumeration  if  we  discover  a  solution  we 
can  prove  is  optimal,  or  acceptably  close  to  optimal  (within  some  moderating  decomposi¬ 
tion  tolerance).  However,  this  algorithm  has  exponential  worst-case  runtime.  If  we  place 
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an  a  priori  limit  on  the  number  of  enumeration  steps,  we  have  a  non-eyelie  heuristie 
seareh  with  domain  limited  by  interdietion  resouree  availability  and  diversity  among  all 
eandidate  delay  plans.  Our  master  problem  objeetive,  whieh  essentially  guides  the 
seareh,  expresses  a  myopie,  greedy  goal:  maximally  delay  the  most-reeently  expedited 
plan.  Other  objeetives  may  provide  faster  termination. 

5,  Finding  Tasks  Affected  by  a  Delay  Plan 

Because  we  assume  the  interdictor’s  actions  to  delay  a  task  will  force  the  prolif- 
erator  to  restart  any  tasks  affected  by  the  delay  plan,  we  identify  those  tasks  before 
proceeding.  A  solution  to  this  model  associates  a  “mark”  with  each  affected  task: 


min  ^f/. 

u  ^ 

(cO) 

i 

s.t.  U.  =  \ 

(cl) 

u,  >  u, 

V(/,7) ^  A- (CHILD SS) 

(c2) 

£/,  E  {0,1} 

Vz  e  A-cr 

(c3) 

The  objective  (cO)  seeks  to  minimize  the  number  of  tasks  marked.  Constraints 
(cl)  mark  all  of  the  tasks  in  the  current  delay  plan  while  each  constraint  (c2)  forces  a  task 
to  be  marked  if  any  of  its  predecessors  are.  Constraints  (c3)  require  the  marks  to  be  bi¬ 
nary.  A  preferable  alternative  to  this  overly  formal  integer  linear  programming 
formulation  is  to  find  the  successor  tasks  using  a  depth-first  search  algorithm  [e.g., 
Abuja,  et  al,  1993,  pp.  73-76].  The  search  algorithm  marks  an  interdicted  task  and  then 
follows  the  path  of  its  successors  until  reaching  the  finish  task,  marking  each  unmarked 
task  visited  along  the  way  and  then  repeating  this  process  for  each  interdicted  task  in  the 
incumbent  delay  plan. 

6.  Expediting  an  Interdicted  Project 

Once  we  have  determined  the  affected  tasks,  we  determine  the  optimal  expediting 
plan  using  updated  task  data. 
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min  S 

S,D,3 

-A 

(dO) 

s.t.  Sj 

+  m(l -£).)> 

V(i,7)  e 

FS 

(dl) 

-  w,)X^, 

-  S,)+ delay  ^x^^yiagy 

1 

IV 

V(i,7)  e 

FF 

(d2) 

s,  +  ((i-a(i 

-  w,)X^, 

-S^+ delay. lag y 

a)K- 

1 

IV 

V(i,7)  e 

SF 

(d3) 

s,  +  lagy 

+  m(l-Z),)> 

Si  +  lagy 

V(i,7)  e 

SS 

(d4) 

s 

r 

E  (1- 

ie^N—rr—r 

Pi{l-U, 

))  (fZ,>A  +  A, )  -  budget  -  budgetUsed 

(d5) 

,)Xff-A 

+  rer^^S.)<a^-d^ 

\/r  e  R 

(d6) 

6, A  <  E 

A  ^ 

bp,  Vz  e  DN 

(d7) 

j\{iJ)GA  -(CHlLDnFF) 

^A 

V(zj)e  A -CHILD 

|ze  (A-OA) 

(d8) 

^A 

e  CHILD  \ 

(d9) 

(3(zz,z)  e  A  1  ((zz,z)  ^ 

CHILD  &(i, 

n)  i  child)) 

p.(not  3(zz,7)  ^  a  \  n 

e  DN) 

J 

0  <  (5^  < 

Vze  N-DN-cj-t 

(dlO) 

A 

e  {0,1} 

yi^N 

(dll) 

A 

=  0,D^=\ 

(dl2) 

Vz  e  A  A  <  t  and  u 

.=0) 

(dl3) 

Vz  e  A  A  <  t  and  u 

,  =  o) 

(dl4) 

>t 

Vz  e  A-5r|M,=l 

(dl5) 

As  in  the  case  with  no  delaying  effort,  the  proliferator  seeks  to  minimize  the  over¬ 


all  length  of  the  project,  expressed  in  (dO).  Constraints  (dl)  through  (d6)  are  similar  to 
constraints  (al)  through  (a6)  from  the  interdiction- free  expediting  model,  but  have  been 
modified  to  account  for  the  effects  of  a  delay  plan  on  the  cost  and  duration  of  each  task. 
The  term  (1  -  /».(!  -  ii^  restores  the  full  cost  and  duration  of  any  task  affected  by  the  de¬ 
lay  plan,  and  accounts  for  the  remaining  cost  and  duration  (the  total  minus  any  cost 
expended  or  progress  achieved  before  time  t)  of  any  task  not  affected  by  the  delay  plan. 
The  budget  and  resource  constraints,  (d5)  and  (d6),  deduct  the  amounts  consumed  up  to 
the  current  time,  computed  with  the  assumption  of  no  interdiction  delays,  from  their  re- 
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spective  budgets.  Constraints  (d7)  through  (dl2)  are  identieal  to  eonstraints  (a7)  through 
(al2). 

Constraints  (al3)  through  (al6)  have  been  removed  and  replaeed  by  (dl3) 
through  (dl5)  that  foree  the  resulting  erashing  plan  to  respeet  the  progress  that  was  made 
before  the  development  projeet  was  diseovered  and  the  delay  plan  formed.  Constraints 
(dl3)  and  (dl4)  fix  the  start  times  and  erashing  effort  of  all  tasks  that  were  in  progress  or 
eomplete  before  the  interdietion  and  are  not  affeeted  by  the  delay  plan.  Constraints  (dl5) 
require  that  all  tasks  affeeted  by  the  delay  plan  start  no  earlier  than  the  eurrent  time.  This 
effeetively  restarts  work  on  any  tasks  that  were  in  progress  before,  but  were  interrupted 
by  the  interdietion.  The  proliferator  is  not  restrieted  from  ehoosing  a  different  deeision 
plan  and  abandoning  (partially)  eompleted  tasks,  if  required. 
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III.  A  COVERT  PROJECT  TO  DEVELOP  A  FISSION  BOMB 


A.  A  CASE  STUDY 

A  rogue  nation  is  suspected  of  pursuing  the  development  of  a  uranium  fission 
weapon.  The  United  States  National  Command  Authority  is  keenly  interested  in  delaying 
this  weapons  program  for  as  long  as  possible.  The  range  of  delaying  tactics  is  consider¬ 
able,  from  diplomatic  pressure,  to  asking  the  United  Nations  to  impose  sanctions,  to 
covert  or  even  overt  military  strikes. 

The  basics  of  nuclear  weapon  design  are  now  well  known  and  publicly  available: 
many  of  the  details  from  the  early  weapons  programs  in  the  United  States  and  elsewhere 
have  been  declassified  and  published  in  the  open  literature.  This  reduces  the  basic  phys¬ 
ics  research  required  to  develop  a  bomb,  but  unless  it  is  possible  to  buy  weapons-grade 
uranium  or  plutonium  on  the  open  market  (not  very  likely  given  the  tight  international 
controls  in  place  today),  developing  an  organic  capability  to  produce  these  materials  re¬ 
quires  a  substantial  investment  in  industrial  infrastructure. 

B.  ASSUMPTIONS 

We  postulate  a  medium-sized  Southwest  Asian  country  with  a  primarily  agrarian 
economy  [Harney,  2003].  The  population  is  generally  well  educated  by  several  modern, 
well-equipped  research  universities.  Most  of  the  nation’s  industrial  capacity  is  concen¬ 
trated  in  a  half-dozen  large  cities.  Many  industries  are  underdeveloped  by  Western 
standards,  but  are  growing  steadily  in  size  and  capability.  In  addition  to  extensive  oil  and 
gas  production,  the  chemical  industries  are  modem  and  well  developed.  The  country  has 
substantial  reserves  of  uranium  ore  and  is  a  well-established  producer  of  both  ore  and 
concentrated  uranium  “yellowcake”  for  the  international  market.  The  country  has  ratified 
the  Nuclear  Nonproliferation  Treaty  and  all  IAEA  safeguards  are  in  place. 

Currently,  50%  of  their  electrical  power  generation  is  from  oil-  or  gas-fired  plants, 
25%  from  hydroelectric  dams  and  25%  from  nuclear  power.  Power  demands  are  growing 
steadily  due  to  population  and  industrial  growth.  Because  the  growth  in  power  demand  is 
geographically  distant  from  the  oil  and  gas  production  centers  and  hydroelectric  produc¬ 
tion  is  already  maximized,  the  most  practical  way  to  satisfy  the  growing  demand  for 
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electricity  is  to  expand  nuclear  generation  capacity.  The  nation  already  has  three  opera¬ 
tional  commercial  nuclear  power  plants  with  two  more  under  construction.  They  also 
have  one  working  research  reactor  affiliated  with  a  university. 

The  military  is  not  large,  but  has  sufficient  capability  to  defend  against  small  re¬ 
gional  threats.  The  army  has  a  number  of  short-  and  medium-range  conventionally  armed 
missiles.  The  air  force  has  no  heavy  bombers,  but  is  equipped  with  a  few  dozen  fighters 
and  medium-range  fighter-bombers.  The  navy  has  a  number  of  small  coastal  defense 
craft.  The  nation  possesses  a  modest  conventional  munitions  production  capability. 

This  country  seeks  its  own  nuclear  weapon  to  counter  growing  threats  from  its 
neighbors.  At  least  three  other  countries  in  the  region  are  known  or  suspected  to  possess 
a  nuclear  arsenal  and  one  nuclear- armed  neighbor  is  growing  more  hostile.  There  is  also 
a  strong  desire  to  be  able  to  use  the  nuclear  weapons  as  a  bargaining  tool  to  obtain  acqui¬ 
escence  and  cooperation  from  several  smaller  neighbors  and  to  gain  credibility  on  the 
international  stage  as  a  key  player  in  the  region.  The  consensus  among  military  and  po¬ 
litical  advisors  is  that  a  total  nuclear  stockpile  of  several  dozen  weapons  will  be  sufficient 
to  achieve  these  goals. 

Because  nuclear  power  generation  is  vital  to  the  country’s  economy  and  grows 
more  so  each  year,  they  have  no  desire  to  immediately  jeopardize  this  by  openly  violating 
the  NPT.  Thus,  any  weapons  program  must  be  covert  and  completely  independent  of  any 
of  the  existing  nuclear  facilities  that  are  being  monitored  by  the  IAEA.  History  shows 
that  most  nations  seeking  nuclear  weapons  have  followed  this  route  and  established  sepa¬ 
rate  military  programs  rather  than  divert  nuclear  material  from  their  safeguarded  civilian 
facilities  [NERAC,  2001]. 

Because  uranium  mining  and  yellowcake  production  are  not  subject  to  NPT  safe¬ 
guards  [U.S.  Congress  OTA,  1993,  p.  137],  the  safest  path  to  a  weapons  program  is 
through  development  of  a  uranium  enrichment  capability  fed  by  existing  yellowcake  pro¬ 
duction.  As  illustrated  in  Eigure  5,  this  will  provide  the  raw  fissile  material  required  for 
either  of  two  types  of  uranium  fission  weapon  (gun  or  implosion)  as  well  as  being  able  to 
produce  fuel  for  the  commercial  reactors  once  their  NPT  violations  have  been  discovered 
and  United  Nations  sanctions  are  imposed.  Given  the  proliferator’s  goal  of  completing 
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Figure  5.  Developing  a  nuclear  weapon  [from  Spears,  2001]. 

We  study  this  up  through  the  weapons  stockpile. 
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his  first  weapon  as  soon  as  possible,  we  assume  that  he  will  first  pursue  a  gun-type 
weapon,  the  same  design  used  by  the  “Little  Boy”  bomb  dropped  on  Hiroshima,  Japan, 
ending  World  War  II.  Little  Boy  was  a  relatively  erude  deviee,  and  its  designers  were  so 
eonfident  in  its  eonstruetion  that  Hiroshima  was  its  first  fiill-seale  test  [Rhodes,  1986  and 
Rhodes,  1995]. 

The  new  nuelear  researeh  and  produetion  faeilities  ean  realistieally  be  hidden  in 
existing  industrial  parks,  where  they  should  eseape  notiee  due  to  the  ongoing  growth  of 
legitimate  industry.  Any  additional  site  seeurity  required  for  the  nuelear  faeilities  ean  be 
hidden  from  outside  observers  in  a  number  of  praetieal  ways. 

C.  DATA  DEVELOPMENT 

To  reduee  the  ehanee  of  deteetion  while  still  aehieving  suffieient  eapaeity  to  meet 
the  arsenal  goal  in  a  reasonable  amount  of  time,  the  faeilities  will  be  designed  to  be  able 
to  produee  five  weapons  a  year.  As  shown  by  the  ealeulations  in  Figures  6  and  7,  this 
produetion  level  will  require  an  annual  input  of  250  kilograms  of  highly  enriehed  ura¬ 
nium  (HEU).  Enriehment  will  thus  require  an  input  of  about  68  metrie  tons  of 
yelloweake  eaeh  year.  Beeause  IAEA  safeguards  for  yelloweake  eover  only  imports  and 
exports,  eovert  diversion  of  this  relatively  small  quantity  (5.6  metrie  tons — one  truek- 
load — a  month)  from  the  existing  produetion  faeilities  should  be  easy. 

Designing  a  nuelear  weapon  and  eonstrueting  the  sophistieated  support  faeilities 
required  to  build  five  weapons  a  year  is  a  daunting  projeet.  Neeessary  aehievements  in- 
elude: 

•  Covert  diversion  of  68  metrie  tons  of  yelloweake  annually; 

•  produetion  of  enriehment  plant  feed  material  (uranium  hexafluoride,  UF() 
from  yelloweake; 

•  uranium  enriehment,  ineluding  the  ehoiee  of  method  to  employ; 

•  eonversion  of  highly  enriehed  UF^,  to  uranium  metal;  and 

•  design  and  eonstruetion  of  the  aetual  weapons. 

Eigure  8  displays  the  tasks  ineluded  in  the  ease  study.  We  assess  the  requirements 
for  speeialized  equipment  from  the  ehemieal  proeesses  shown  in  Eigure  6.  There  are  a 
total  of  194  tasks  (nodes)  and  583  task  partial  orders  (ares)  in  the  resulting  projeet  net¬ 
work.  The  proliferator  must  manage  five  non-renewable  resourees:  energy,  materials, 
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Technology 

Separation 

Factor 

No,  Of  Stages 
For  90%  HEU 

KWh/SWU 

KWFor 
50,000  SWUA^r 

Gaseous  diffusion 

1.004-1.0045 

3,500-4,000 

2,500 

14,270 

Gas  eentrifuge 

1.2-1.5 

40-90 

100-200 

571-1,142 

Aerodynamie 

1.015-1.030 

540-1100 

3,600-4,000 

20,548-22,831 

Table  1.  Efficiency  of  uranium  enrichment  technologies  [after  U,S,  Congress  OTA, 
1993,  p.  143]. 

Regardless  of  the  teehnology  used,  produeing  HEU  from  natural  uranium  re¬ 
quires  multiple  equipment  stages  arranged  in  a  progressive  enriehment 
easeade.  The  separation  factor  is  the  ratio  of  the  relative  enriehment  (17235  to 
17238)  of  the  eoneentrated  produet  to  that  of  the  depleted  tails  from  the  output 
of  any  one  stage  of  the  easeade.  The  number  of  stages  required  to  produee 
HEU  assumes  the  final  tails  eontain  less  than  0.3%  17235.  Kilowatt-hours  per 
SWU  measures  the  amount  of  energy  required  to  produee  approximately  5 
grams  of  HEU.  The  final  eolumn  is  the  average  power  eonsumption  rate  of  an 
enriehment  easeade  produeing  about  250  kg  of  HEU  a  year,  enough  to  eon- 
struet  five  fission  weapons. 

professional  labor  (e.g.,  seientists  and  engineers),  skilled  labor  (e.g.,  maehinists)  and  un¬ 
skilled  labor.  He  must  ehoose  one  and  only  one  enriehment  teehnology  to  pursue.  His 
ehoiees  are  gas  eentrifuges,  gaseous  diffusion  or  aerodynamie  enriehment.  Table  1  lists 
information  used  to  determine  the  required  number  of  enriehment  maehines  and  energy 
eonsumption  for  eaeh  teehnology.  A  separative  work  unit  (SWU)  is  a  measure  of  the  ef¬ 
fort  required  to  separate  17235  and  17238  isotopes  during  enriehment.  Produetion  of  1  kg  of 
HEU  from  natural  uranium  requires  approximately  200  SWU  [Nonproliferation  Poliey 
Edueation  Center,  2004]. 

Our  ease  study  employs  only  mild  interdietions,  the  better  to  provide  a  nuaneed 
solution  exereise.  It  is  not  lost  on  us  that  a  single  Tomahawk  missile  ean  render  an  en¬ 
riehment  eentrifuge  gallery  inoperable,  and  essentially  unoeeupiable  for  millennia. 
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Feed  stock  preparation  (yellowcake  to  uranium  hexafluoride): 

Stainless  steel  vessel  (dissolution  of  yelloweake  in  nitrie  aeid): 

+  6  HNO^  3  UO^ {NO, )^+2H,0  +  H, 

842 +  6x  63  =  3x394 +  2x18  +  2 
Stainless  steel  boiler  (boil  down  of  nitrate  solution): 

UO,{NO,)^  +  6  H,0  >  UO,{NO,)^  * 6  H,0 

394  +  6x18  =  502 

High- temperature  stainless  steel  boiler  (thermal  deeomposition  of  nitrate): 

UO,{NO, * 6  H,0  > U0,+6H,0  +  2  NO^ 

502  =  286  +  6x18  +  108 
Gas-solid  reaetor  vessel  (reduetion  of  uranium  trioxide): 

UO,  +  H,  >  UO,  +  H,0  or  3UO,+2  + N,+3  H,0 

286  +  2  =  270  +  18  or  3x286  +  2x17  =  3x270  +  28  +  3x18 

Stainless  steel  reaetion  vessel: 

UO,  +4HF - >UF,  +2H,0 

270  +  4x20  =  314  +  2x18 

Ultrahigh-temperature  gas-solid  reaetor  vessel  (produetion  of  uranium  hexafluoride): 

UF,  +  F,  ^ 

314  +  38  =  352 

Material  conversion  (uranium  hexafluoride  to  metal): 

Gas-phase  reaetor  with  partieulate  separation: 

GFg  +  H,  UF^  (5)  +  2HF 

352  +  2  =  314  +  2x20 
High-temperature  metallurgieal  furnaee: 

UF^  (5)  +  2  Ca  U{l)  +  2  CaF,  or  UF,  (5)  +  2  Mg  u{l)  +  2  MgF, 

_ 314  +  2x40  =  238  +  2x78  or  314  +  2x24.3  =  238  +  2x62.3 _ 

Figure  6,  Uranium  enrichment  chemistry  [from  Harney,  2003]. 

We  study  these  ehemieal  proeesses  to  deduee  the  faeilities  and  equipment  re¬ 
quired  to  produee  enriehed  uranium  metal  from  yelloweake.  The  numbers 
beneath  the  ehemieal  reaetion  formulas  are  moleeular  weights. 
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1  kg  HEU  metal  requires  1 .302  kg  HEUF4  and  (0.336  kg  Ca  or  0.204  kg  Mg) 

1  kg  HEUF4  requires  1.122  kg  HEUF(,  and  0.0064  kg  H2 
1 .302  kg  HEUF4  requires  1 .461  kg  HEUF(,  and  0.0084  kg  H2 

1  kg  7/2  can  be  produeed  from  5.667  kg  of  NFIt, 

1  kg  HEUF(,  requires  232  kg  UF(, 

1.461  kg  HEUF(,  requires  339  kg  UF(, 

1  kg  UF(,  requires  0.892  kg  UF4  and  0.108  kg  F2 
339  kg  UF(,  requires  302  kg  UF4  and  36.6  kg  F2 

1  kg  UF4  requires  0.860  kg  UO2  and  0.255  kg  FIE 
302  kg  UF4  requires  260  kg  UO2  and  77.0  kg  HF 

1  kg  UO2  requires  1.059  kg  UO2  and  0.0074  kg  7/2 
260  kg  UO2  requires  275  kg  UO2,  and  1.924  kg  FI2 

1  kg  UO2,  requires  1.755  kg  of  hexahydrate 
275  kg  UO2,  requires  483  kg  of  hexahydrate 

1  kg  of  hexahydrate  requires  0.559  kg  yelloweake  and  0.251  kg  HNOj, 

483  kg  of  hexahydrate  requires  270  kg  yelloweake  and  121  kg  HNOj, 

Eaeh  weapon  requires  50  kg  93%  HEEl  metal 

5  weapons/yr  requires  250  kg  HEEl  metal 

250  kg  HEU  requires;  67,500  kg  yelloweake 

30.250  kg  HNO3 

2.1  kg  Hi  (or  11.9  kg  NHi) 

19.250  kgEfF 
9,150  kg  Fi 

84  kg  Ca  or  51  kg  Mg 

Figure  7.  Calculation  of  annual  raw  material  input  requirements  [after  Harney, 
2003]. 

We  assess  the  quantity  of  raw  materials  required  to  sustain  produetion  of  five 
weapons  per  year.  Requirements  are  eomputed  from  the  ehemieal  reaetions 
and  moleeular  weights  shown  in  Eigure  6. 
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•  Nuclear  Weapon  Development 

■  Divert  Commercial  Yellowcake 

o  Design  yellowcake  plant  modifications 
o  Modify  yellowcake  plant 
o  Divert  yellowcake 

■  Produce  Enrichment  Plant  Feed  Material 

o  Design  fluoridation  plant 
o  Acquire  plant  site 
o  Prepare  plant  site 
o  Acquisition  Of  Plant  Components 

•  Stainless  steel  mixing  vessel 

•  Distilled  water  system 

•  Nitric  acid  storage  tank 

•  Stainless  steel  boiler 

•  Thermal  decomposition  vessel 

•  Drying  kiln 

•  Gas-solid  high-temperature  reac¬ 
tion  vessel 

•  Hydrogen  gas  (or  ammonia)  stor¬ 
age  tank 

•  Stainless  steel  reaction  vessel 

•  Hydrogen  fluoride  storage  tank 

•  Gas-solid  ultra  high-temperature 
reaction  vessel 

•  Fluorine  storage  tank 

•  Hexafluoride  condensing  vessel 

•  Pumps  and  piping 

o  Assembly  and  integration 
o  Operation 

■  Uranium  Enrichment 

o  Choose  enrichment  technology 
o  Gas  Centrifuge 

•  Design  basic  gas  centrifuge 

•  Acquire  Research  Components 

■  Rotor  tubes 

■  Air  bearing  systems 

■  Motors 

■  End  caps 

■  Centrifuge  cases 

■  Pumps  and  piping 

•  Assemble  research  centrifuges 

•  Test  and  evaluation  of  research 
centrifuges 

•  Design  production  centrifuges 

•  Design  enrichment  cascade 

•  Design  enrichment  plant 

•  Acquire  plant  site 

•  Prepare  plant  site 

•  Acquire  Production  Components 

■  Rotor  tubes 

■  Air  bearing  systems 

■  Motors 

■  End  caps 

■  Centrifuge  cases 

■  Pumps  and  piping 

•  Assemble  production  centrifuges 

•  Integrate  enrichment  cascade 
o  Gas  Diffusion 


•  Design  basic  gas  diffusion  system 

•  Acquire  Research  Components 

■  Diffusion  barriers 

■  Heat  exchangers 

■  Pumps  and  piping 

•  Assemble  research  devices 

•  Test  and  evaluation  of  research  de¬ 
vices 

•  Design  production  gas  diffusion 
devices 

•  Design  enrichment  cascade 

•  Design  enrichment  plant 

•  Acquire  plant  site 

•  Prepare  plant  site 

•  Acquire  Production  Components 

■  Diffusion  barriers 

■  Heat  exchangers 

■  Pumps  and  piping 

•  Assemble  production  devices 

•  Integrate  enrichment  cascade 
o  Aerodynamic 

•  Design  basic  aerodynamic  enrich¬ 
ment  device 

•  Acquire  Research  Components 

■  Vortex  unit 

■  Pumps  and  piping 

•  Assemble  research  devices 

•  Test  and  evaluation  of  research  de¬ 
vices 

•  Design  production  devices 

•  Design  enrichment  cascade 

•  Design  enrichment  plant 

•  Acquire  plant  site 

•  Prepare  plant  site 

•  Acquire  Production  Components 

■  Vortex  unit 

■  Pumps  and  piping 

•  Assemble  production  devices 

•  Integrate  enrichment  cascade 
o  Cascade  loading 

o  Produce  enriched  and  depleted  material 
■  Prepare  Uranium  Metai 
o  Design  metal  plant 
o  Acquire  plant  site 
o  Prepare  plant  site 

o  Acquire  Piant  Components  (Enriched 
Metai) 

•  Gas-phase  reactor  with  particulate 
collection 

•  Hydrogen  storage  tank 

•  Metallurgical  furnace 

•  Hafnia  crucibles 

o  Acquire  Piant  Components  (Depieted 
Metai) 

•  Gas-phase  reactor  with  particulate 
collection 

•  Hydrogen  storage  tank 

•  Metallurgical  furnace 
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•  Hafnia  crucibles 

•  Double-base  propellant  powder 

o 

Integrate  components 

•  Polonium 

o 

Produce  natural  uranium  metal 

•  Beryllium  powder 

o 

Produce  depleted  and  enriched  uranium 

•  Detonator  and  explosive  train  com- 

metal 

ponents 

■  Prepare  Nuclear  Explosive  Devices 

o  Fabricate  Research  Device  Compo- 

o 

Design  Gun  Device  Components 

nents  (Enriched  Uranium  Prototype) 

•  Gun 

•  Gun  barrel 

•  Propellant 

•  Breech  mechanism 

•  Critical  Core 

•  Cast  enriched  uranium  components 

■  Fissionable  receiver 

•  Cast  depleted  uranium  tamper 

■  Fissionable  projectile 

•  Machine  enriched  uranium  receiver 

■  Tamper 

•  Machine  enriched  uranium  projec- 

■  Initiator 

tile 

•  Safety  and  arming  devices 

•  Machine  depleted  uranium  tamper 

•  Fuse 

•  Initiator 

o 

Design  weapon  assembly  plant 

•  Propellant  charge 

o 

Acquire  plant  site 

•  Detonator  and  explosive  train 

o 

Prepare  plant  site 

o  Assemble  research  devices  (enriched 

o 

Acquire  Fabrication  Devices 

uranium  prototype) 

•  Large -diameter  precision  lathe 

o  Test  Research  Devices  (Sub-Critical) 

•  Inert-gas  environment  precision 

•  Verify  critical  mass 

milling  machine 

•  Verify  gun  velocity 

•  Metallurgical  furnace 

•  Delivery  vehicle  compatibility 

•  Hafnium  cracibles 

mock-up 

•  Inert-gas  environment  casting  sys- 

o  Test  full-scale  device 

tern 

o  Finalize  production  weapon  design 

o 

Acquire  Research  Device  Compo- 

o  Acquire  Weapon  Components 

nents  (Natural  Uranium  Prototype) 

•  High-strength  steel  cylinder 

•  High-strength  steel  cylinder 

•  Double-base  propellant  powder 

•  Double-base  propellant  powder 

•  Polonium 

•  Polonium 

•  Beryllium  powder 

•  Beryllium  powder 

•  Detonator  and  explosive  train  com- 

•  Detonator  and  explosive  train  com- 

ponents 

ponents 

o  Fabricate  Weapon  Components 

o 

Fabricate  Research  Device  Compo- 

•  Gun  barrel 

nents  (Natural  Uranium  Prototype) 

•  Breech  mechanism 

•  Gun  barrel 

•  Cast  enriched  uranium  components 

•  Breech  mechanism 

•  Cast  depleted  uranium  tamper 

•  Cast  uranium  components 

•  Machine  enriched  uranium  receiver 

•  Cast  uranium  tamper 

•  Machine  enriched  uranium  projec- 

•  Machine  uranium  receiver 

tile 

•  Machine  uranium  projectile 

•  Machine  depleted  uranium  tamper 

•  Machine  uranium  tamper 

•  Initiator 

•  Initiator 

•  Propellant  charge 

•  Propellant  charge 

•  Detonator  and  explosive  train 

•  Detonator  and  explosive  train 

•  Fuse 

o 

Assemble  research  devices  (natural  ura- 

•  Safety  and  arming  device 

nium  prototype) 

•  Weapon  case  and  structure 

o 

Acquire  Research  Device  Compo- 

o  Assemble  weapon  components 

nents  (Enriched  Uranium  Prototype) 

•  High-strength  steel  cylinder 

o  Production  weapons  deliveries 

Figure  8.  Fission  weapon  development  tasks  used  in  our  case  study  [after  Harney, 
2003]. 

Summary  tasks  are  indicated  in  bold  and  the  Choose  enrichment  technology 
decision  node  is  in  italics. 
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D,  INSIGHTS 

We  present  our  enumeration  proeedure  as  a  tool  that  provides  a  quantitative  eom- 
parison  of  the  myriad  options  available  to  an  interdietor  for  senior  nonproliferation  poliey 
makers.  Analysis  of  an  interdietor’s  options  in  our  ease  study  ean  take  several  direetions. 
Ignoring  for  the  moment  our  two-sided  optimization,  the  most  straightforward  applieation 
answers  “what-if  ’  questions  manipulating  just  the  proliferator’s  planning  problem.  Pos¬ 
tulate  a  delay  plan  based  on  eurrent  or  proposed  national  poliey,  infliet  it  on  the  projeet 
and  immediately  visualize  and  evaluate  the  proliferator’s  optimal  expediting  plan  in  re¬ 
sponse  to  that  delay  plan.  Planners  ean  vary  input  parameters  sueh  as  the  expediting 
budget,  resouree  availabilities,  et  cetera,  over  ranges  of  uneertainty  to  test  the  sensitivity 
of  the  expediting  plan.  They  might  also  alter  the  plan  to  enhanee  its  eonsisteney  with 
telltale  symptoms  of  its  aetivities,  to  better  assess  the  true  stage  of  proliferation.  This  is 
useful,  but  gives  no  information  about  whether  a  better  delay  plan  exists. 

Using  our  two-sided  enumeration  to  find  the  best  delay  plan  for  a  single  eombina- 
tion  of  input  parameters,  we  ean  reeover  the  expediting  plan  resulting  from  any  eandidate 
delay  plan  that  is  diseovered.  This  ease  study  reveals  many  instanees  of  multiple  delay 
plans  that  all  produee  the  same  maximum  expedited  projeet  length.  A  set  of  sueh  delay 
plans  ean  be  presented  to  poliey  makers  for  a  more  subjeetive  evaluation. 

The  final  mode  of  study  is  also  the  most  powerful.  We  again  employ  our  enu¬ 
meration  proeedure  to  find  delay  plans,  but  we  now  allow  input  parameters  to  vary  over 
ranges  of  uneertainty.  The  results  probe  the  sensitivity  of  the  ehoiee  of  delay  plan  to  the 
inputs  as  well  as  the  proliferator’s  response  to  that  delay  plan.  The  relative  frequency 
with  which  a  particular  task  appears  in  the  optimal  delay  plans  turns  out  to  be  a  key  in¬ 
dicator  of  vulnerability  and  robustness.  If  one  partieular  task  appears  in  nearly  every 
optimal  delay  plan,  regardless  what  the  input,  that  is  a  elear  signal  of  vulnerability. 

Our  proliferator  has  an  essentially  unlimited  supply  of  yelloweake  at  his  disposal; 
therefore,  our  ease  study  does  not  inelude  the  possibility  of  interdieting  the  uranium  pro- 
duetion.  When  we  allow  unlimited  quantities  of  the  other  resourees,  our  proliferator 
ehooses  aerodynamie  enriehment  and  eompletes  his  projeet  in  260  weeks  (five  years), 
given  that  we  do  nothing  to  delay  him.  If  we  assume  the  proliferator  has  no  knowledge 
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—  340  weeks:  Covert  interdiction 

16  weeks 

- 

—  324  weeks:  Complete  transparency 

64  weeks  ^ 

-  260  weeks:  Covert  proliferation,  no  interdiction 

Figure  9.  Value  of  information  to  the  proliferator  and  the  interdictor. 

The  proliferator  achieves  his  goal  64  weeks  (over  a  year)  sooner  by  conceal¬ 
ing  his  project  from  the  interdictor.  The  interdictor  extends  the  project  16 
weeks  when  the  interdiction  is  covert. 

of  our  intention  to  interdict  his  project,  to  maximize  the  total  delay  we  choose  tasks  to 
interdict  that  are  on  his  critical  path.  By  choosing  the  two  tasks  on  the  critical  path  with 
the  greatest  delay,  we  can  extend  the  proliferator’ s  project  to  340  weeks  (6X  years,  a 
31%  increase),  assuming  he  does  not  discover  or  respond  to  our  interdiction.  To  achieve 
this  delay,  we  must  interdict  cascade  loading  and  acquisition  of  either  the  vortex  unit  or 
pumps  and  piping  for  the  research  enrichment  machines. 

With  an  intelligent  adversary,  our  two-sided  shows  that  even  if  the  proliferator  is 
fully  aware  of  our  intention  to  interdict  his  project  and  which  tasks  we  will  delay,  the  best 
he  can  do  is  complete  the  project  in  324.4  weeks  (about  6X  years).  This  is  a  25%  in¬ 
crease  over  the  length  of  the  project  when  we  do  nothing.  Even  though  the  proliferator 
knows  the  project  will  be  interdicted  and  has  free  reign  to  compensate,  he  is  only  able  to 
save  15.6  weeks  or  4.6%  of  the  project  length  we  achieve  by  interdicting  with  complete 
surprise.  The  two  tasks  in  the  optimal  delay  plan  are  cascade  loading  and  acquisition  of 
pumps  and  piping  for  the  production  aerodynamic  enrichment  machinery.  However,  be¬ 
cause  we  interdict  equipment  required  for  aerodynamic  enrichment,  the  proliferator 
chooses  to  save  time  by  using  gaseous  diffusion  instead. 

Figure  9  shows  the  value  of  secrecy.  By  maintaining  a  covert  project  that  cannot 
be  interdicted  because  no  one  knows  it  exists,  the  proliferator  can  complete  his  first  nu¬ 
clear  weapon  almost  15  months  sooner  than  if  he  anticipates  being  interdicted.  The 
interdictor  can  delay  the  project  16  weeks  longer  if  he  achieves  surprise. 
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Cascade  loading  appears  in  all  of  these  delay  plans.  In  every  ease,  with  both 
eomplete  surprise  and  eomplete  transpareney,  it  is  the  most  vulnerable  task  in  this  pro- 
jeet.  This  is  useful  information  to  a  poliey  maker  or  military  planner  eontemplating  the 
obstruetion  of  the  proliferator’s  projeet. 

Figure  10  illustrates  the  effeet  on  the  expedited  projeet  length  of  varying  only  one 
input  parameter.  There  is  a  elear  ehange  in  the  projeet  length  when  the  materials  budget 
is  somewhere  between  $34  and  $35  million.  Inspeetion  of  the  resulting  expediting  plans 
reveals  that  the  reason  for  the  step  is  a  ehange  in  the  deeision  plan.  With  a  small  materi¬ 
als  budget,  the  proliferator  is  foreed  to  ehoose  aerodynamie  enriehment,  extending  the 
length  of  the  projeet.  Onee  a  suffieient  materials  budget  is  provided,  he  ean  ehoose  gas- 
diffusion  in  order  to  finish  his  first  weapon  more  quiekly.  Figure  1 1  magnifies  the  transi¬ 
tion  in  Figure  10.  Infieetion  points  like  this  are  easy  to  diseover  by  systematioally 
searehing  the  domain  of  resouree  availability. 

Figure  12  shows  the  effeet  of  varying  both  the  expediting  budget  and  the  delaying 
budget  simultaneously.  It  is  elear  that  no  matter  what  the  interdietor  does,  inereasing  the 
expediting  budget  available  to  the  proliferator  reduees  the  time  needed  to  finish  the  pro¬ 
jeet.  Onee  the  delay  budget  reaehes  a  eertain  point,  there  are  no  more  expensive  delaying 
aetions  available  to  the  interdietor  that  eause  more  delay  than  the  inexpensive  one  identi¬ 
fied  with  the  lower  interdietion  budget. 

Figure  13  shows  the  effeet  of  varying  the  expediting  budget  and  the  number  of  in- 
terdietions  allowed  without  eonsidering  the  interdietion  budget.  Logieally,  allowing 
more  interdietions  admits  more  delay  to  the  projeet. 
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Project  Length  versus  Materials  Budget 


Materials  Budget  ($M) 

Figure  10,  Project  length  versus  materials  budget. 

The  proliferator  is  allowed  a  $100M  crashing  budget,  400  man-months  of  pro¬ 
fessional  labor  and  unrestricted  amounts  of  the  other  resources.  The 
interdictor  is  allowed  two  interdictions  and  $1M  delay  budget.  The  step  be¬ 
tween  $34  and  $35  million  reflects  a  change  in  the  decision  plan.  Lower 
materials  budgets  force  the  expeditor  to  use  aerodynamic  enrichment  while 
higher  budgets  encourage  use  of  gas-diffusion  to  reduce  the  overall  length  of 
the  project  by  3.4  weeks. 

Project  Length  versus  Materials  Budget 


Figure  11,  Project  length  versus  materials  budget,  a  closer  look, 

A  closer  look  at  the  decision  plan  transition  in  Figure  10  shows  that  the 
change  actually  occurs  when  the  materials  budget  is  between  $34.6  and  $34.7 
million. 
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Project  Length  versus  Expediting  Budget  and  Delay  Budget 


Expediting  Budget  ($M)  1 10 


Figure  12,  Project  length  versus  budget  with  one  interdiction. 

Allowing  more  expensive  interdiction  actions  does  not  always  cause  more  de¬ 
lay  to  the  project.  Once  the  delay  budget  reaches  $500K,  no  further 
improvement  in  the  project  delay  is  achieved.  No  resource  is  exhausted  here. 


Project  Length  versus 

Expediting  Budget  and  Number  of  Interdictions 


Number  of 
Interdictions 


Figure  13,  Project  length  versus  expediting  budget  and  number  of  interdictions. 

No  resource  restriction  is  imposed. 
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IV.  IMPLEMENTATION 


We  want  to  show  how  to  solve  praetieal  problems,  so  we  have  integrated  our  in- 
terdietion  proeedure  with  an  off-the-shelf  projeet  management  tool  to  manage  and  display 
the  data.  The  interdietion  optimization  has  been  implemented  using  the  General  Alge- 
braie  Modeling  System  (GAMS)  [Brooke,  et  al,  1998]  and  XA  solver  from  Sunset 
Software  Teehnology  [2003].  The  data  is  stored  and  manipulated  in  a  Mierosoft  Projeet 
(MS  Projeet)  [Mierosoft  Corporation,  2003a]  database.  A  eustom  Visual  Basie  for  Ap- 
plieations  (VBA)  [Mierosoft  Corporation,  2003b]  user  interfaee  hosts  the  optimization 
model  in  GAMS.  The  VBA  eode  and  eustom  display  formatting  are  saved  into  an  MS 
Projeet  template  file. 

By  using  off-the-shelf  software  like  MS  Projeet,  we  gain  an  inexpensive,  flexible 
and  easy  to  use  database  and  display  tool.  As  part  of  the  Mierosoft  Offiee  [Mierosoft 
Corporation,  2003a]  suite  of  business  software,  it  is  widely  available  and  already  familiar 
to  a  large  user  base.  However,  MS  Projeet  does  not  have  the  sophistieated  optimization 
eapabilities  that  we  need.  Like  all  Mierosoft  Offiee  applieations,  it  permits  extensive 
eustomization  with  user-defined  graphieal  user  interfaee  elements  and  VBA  eode.  This 
feature  allows  us  to  easily  hide  the  powerful  optimization  tools  from  the  user  and  instead 
present  only  a  straightforward  graphieal  user  interfaee  to  eontrol  the  optimization  proee- 
dure  and  display  the  results. 

This  ehoiee  of  software  tools  does  present  some  ehallenges,  however.  The  most 
signifieant  of  these  is  the  way  MS  Projeet  requires  a  hierarehieal  strueture  for  the  tasks  in 
the  projeet  to  ease  the  organization  of  large  projeets.  MS  Projeet  displays  eaeh  hierarehi¬ 
eal  task  in  a  Gantt  ehart  [Gantt,  1903]  in  the  style  shown  in  Figure  14.  The  blaek  bar 
represents  the  top-level  summary  task  (the  parent)  and  the  subordinate  tasks  are  its  ehil- 
dren.  Summary  tasks  ean  be  given  their  own  preeedenee  relationship(s)  with  other  tasks 
in  the  projeet.  Aeeommodating  them  requires  introdueing  dummy  preeedenee  ares  into 
the  network  to  ensure  a  summary  task  starts  no  later  than  any  of  its  ehildren  and  to  pre- 
elude  finishing  until  all  of  its  ehildren  finish.  (In  our  model,  this  requires  extra 
eonstraints  for  the  binary  deeision  variables.)  The  summary  task  must  be  eompleted  if 
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any  of  its  children  are  eompleted  and  the  ehildren  must  be  eompleted  if  the  summary  task 
is  forced  to  start  by  other  preeedence  eonstraints. 


S  1  F 

S  2  F 

Figure  14,  Gantt  chart  representation  of  a  hierarchical  summary  task. 

The  summary  task  begins  as  soon  as  any  subordinate  task  starts  and  finishes 
when  all  subordinates  are  finished. 

Figure  15  illustrates  the  eustomized  data  sheet  that  is  used  to  enter  most  of  the 
projeet  data  within  MS  Projeet.  The  hierarehieal  organization  of  the  projeet  tasks  is  evi¬ 
dent  in  this  projeet  view.  This  sheet  also  displays  the  task  durations  and  the  deeision  and 
delay  plans  determined  by  the  interdietion  proeedure.  MS  Projeet  ealeulates  the  eost  of 
eaeh  task  based  on  the  task  duration  and  the  resource  requirements  and  eosts  in  the  pro¬ 
jeet  database.  The  data  sheet  is  adjaeent  to  the  Gantt  ehart  (Figure  16),  enabling  quiek 
switehing  between  the  tabular  data  and  a  graphieal  representation  of  the  task  relationships 
within  the  projeet.  Figure  17  shows  how  tasks  in  the  delay  plan  are  represented  on  the 
Gantt  ehart. 

The  eustom  user  interfaee  is  implemented  with  VBA  eode  distributed  among  two 
user  forms  (eustom  dialog  boxes)  and  a  eode  module.  The  first  user  form  (Figure  18) 
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Duration 

Interdiction 
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Interdiction  Cost 

Decision 

Completed 

Interdicted 
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Cascade  loading 

0% 

122.97 

wks 

$8,1 83,302.11 

112wks 

56  wks 

56  wks 
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Production  of  enriched  and  depleted  material 

0% 
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$300,000.00 

24  wks 

12  wks 

12  wks 

$1 .700,000.00 
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B  Preparation  of  uranium  metal 

0% 

300.94 

wks 

$5,474,500.00 

0  wks 

0  wks 

0  wks 

$0.00 
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Design  of  metal  plant 

0% 

56  wks 

$2,100,000.00 

56  wks 

28  wks 

28  wks 

$450,000.00 
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Acquisition  of  plant  site 

0% 

24  wks 

$140,000.00 

24  wks 

12  wks 

12  wks 

$201 ,000.00 

95 

Preparation  of  plant  site 

0% 
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0% 
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$188,000.00 

0  wks 

0  wks 

0  wks 

$0.00 

97 

Gas-phase  reactor  with  particulate  collection 

0% 

48  wks 

$125,000.00 
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Hydrogen  Storage  Tank 

0% 

48  wks 

$500.00 

48  wks 

24  wks 

24  wks 

$201 ,000.00 
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99 

Metallurgical  Furnace 

0% 

48  wks 

$50,000.00 

48  wks 

24  wks 

24  wks 

$201 ,000.00 
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100 

Hafnia  Crucibles 

0% 

48  wks 

$12,500.00 

48  wks 

24  wks 

24  wks 

$201 ,000.00 

a 

Figure  15,  Selected  project  data  from  our  case  study. 


The  “Deeision”  eolurnn  identifies  deeision  nodes  with  a  checkmark.  The 
“Completed”  eolurnn  indieates  the  deeision  plan  and  the  “Interdieted”  eolurnn 
identifies  the  delay  plan  that  generated  the  displayed  projeet  sehedule. 


34 


Figure  16,  Gantt  chart. 

Vertical  lines  delimit  four-week  time  blocks.  Each  bar  with  the  diagonal  slash 
pattern  indicates  when  a  task  would  be  conducted  if  there  is  no  delay  effort. 
The  shaded  box  above  it  shows  the  effect  of  an  interdiction  delay  elsewhere  in 
the  project.  The  task  bars  extended  with  the  vertical  shading  highlight  time 
saved  by  crashing  those  tasks.  For  example,  the  uranium  component  casting 
task  has  been  delayed  by  58  weeks  and  extra  resources  have  been  committed 
to  reduce  its  duration  by  four  weeks. 

collects  parameters  such  as  the  number  of  iterations  desired  and  the  budget,  creates  data 
files  to  be  included  in  the  GAMS  program  and  then  calls  GAMS.  Each  input  parameter, 
including  the  individual  resource  budgets  (Figure  19),  may  be  given  as  a  range  of  values 
and  the  GAMS  program  will  automatically  solve  the  problem  multiple  times  and  generate 
a  comma-separated  values  (“CSV”)  file  that  can  be  analyzed  using  a  spreadsheet  or 
mathematical  charting  and  analysis  software.  The  second  user  form  (Figure  20)  opens 
automatically  when  GAMS  terminates  or  the  user  presses  the  “Read  Results”  button  in 
the  first  user  form.  It  displays  a  list  containing  the  length  of  the  project  computed  in  each 
iteration  for  each  combination  of  input  parameters.  Using  a  quick  sort  subroutine  adapted 
from  VBA-programmer.com  [2004],  the  user  can  sort  the  list  by  the  values  in  any  col¬ 
umn.  Any  intermediate  result  may  be  selected  and  read  back  into  the  MS  Project 
database  for  display.  The  code  module  provides  links  to  operating  system  services  [Mi¬ 
crosoft  Corporation,  2004]  that  permit  the  optimization  software  to  be  called 
synchronously,  allowing  the  results  to  be  displayed  as  soon  as  GAMS  terminates. 
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‘Production  of  enriched  and  depleted  material 


Figure  17,  An  interdicted  task  on  the  Gantt  chart. 

The  lightly  shaded  section  of  the  task  bar  between  the  dark  shading  and  the 
vertical  shading  shows  the  delay  imposed  on  this  task  by  the  interdiction  ef¬ 
fort.  The  cascade  loading  task  has  thus  been  extended  by  56  weeks  because  of 
the  interdiction. 


Each  instance  of  this  case  study  requires  less  than  one  second  to  generate,  opti¬ 
mize  and  return  to  populate  results  in  MS  Project.  The  subproblem  uses  193  binary  and 
385  continuous  variables  and  1,095  constraints.  XA  solves  this  problem  to  optimality  in 
less  than  0.25  second  on  a  Pentium  4  based  desktop  personal  computer.  The  master  prob¬ 
lem  uses  163  binary  variables  and  two  constraints  plus  one  diversity  constraint  for  each 
prior  delay  plan.  Optimal  solution  times  are  similar  to  the  subproblem. 


Figure  18,  Control  panel  used  to  enter  parameter  ranges  for  the  optimization. 

Clicking  the  “Restrict  solution...”  check  box  runs  the  subproblem  to  evaluate 
a  particular  delay  plan  rather  than  finding  an  optimal  plan.  This  dialog  box 
saves  all  user  input  and  automatically  populates  the  input  fields  when  opened. 


o 

Resource  Name 

Material  Label 

Std.  Rate 

Min  Units 
Available 

Max  Units 
Available 

Units 

Available 

Increment 

1 

Materials 

$K 

$1 ,000.00 

34000 

35000 

100 

2 

Energy 

MWhr 

$10.00 

100000 

100000 

1 

3 

Professional  Labor 

Mmo 

$60,000.00 

1000 

1000 

1 

4 

Skilled  Labor 

Mmo 

$24,000.00 

10000 

10000 

1 

5 

Unskilled  Labor 

Mmo 

$6,000.00 

1000 

1000 

1 

Figure  19,  Resource  data. 

The  “Resource  Sheet”  view  in  MS  Project  is  used  to  enter  the  ranges  for  each 
individual  resource  budget. 
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Figure  20,  Dialog  box  used  to  load  GAMS  results  into  MS  Project, 

The  top  list  displays  every  combination  of  input  parameters  evaluated.  The 


bottom  list  displays  results  of  evaluating  every  candidate  delay  plan  for  the  pa¬ 
rameter  combination  selected  in  the  top  list.  This  example  illustrates  that 
many  different  delay  plans  may  evaluate  to  the  same  overall  project  length. 
The  header  buttons  choose  which  column  to  sort  and  the  arrow  button  changes 
the  sort  direction. 


This  case  study,  along  with  the  MS  Project 
are  available  from  the  author  or  thesis  advisor  upon 


template  file  and  GAMS  source  code 
substantiated  request. 
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V.  WHAT  WE  HAVE  ACCOMPLISHED 


Building  on  decades  of  project  management  research,  we  present  a  novel  and 
flexible  model  that  captures  many  of  the  issues  involved  in  expediting  a  sophisticated  re¬ 
search  and  development  project.  Using  this  tool  as  a  forecast  of  proliferation  plans,  we 
next  seek  an  optimal  set  of  tasks  in  the  project  to  interdict  in  order  to  inflict  the  most  de¬ 
lay  in  its  completion. 

Our  prototypic  decision  support  system  is  implemented  using  off-the-shelf  soft¬ 
ware.  All  user  interaction  and  data  management  are  done  through  software  already 
familiar  to  millions.  Computational  experience  with  our  case  study  shows  that  good  solu¬ 
tions  can  be  found  quickly  (at  less  than  one  second  per  iteration  on  a  desktop  personal 
computer).  “What-if’  questions  regarding  particular  delay  plans  can  be  answered  in  a 
fraction  of  a  second.  An  analyst  can  take  advantage  of  this  speed  to  explore  a  wide  range 
of  assumptions  and  find  the  tasks  that  are  most  critical  to  the  expeditious  completion  of 
the  proliferator’s  project. 

This  low-cost  tool  can  be  used  today.  It  provides  quantitative  information  that 
can  be  used  to  guide  policy  makers  working  to  stop  the  proliferation  of  nuclear,  chemical, 
biological  and/or  radiological  dispersion  weapons,  or  to  disrupt  more  conventional  strate¬ 
gic  industrial  and  commercial  activities. 

Many  research  opportunities  in  project  network  interdiction  remain.  Modifica¬ 
tions  to  the  master  problem  objective  function  may  yield  a  faster  search  for  promising 
candidate  delay  plans.  The  design  of  the  enumeration  procedure  allows  a  tremendous 
amount  of  flexibility  with  the  subproblem  and  the  master  problem.  We  are  limited  only 
by  what  we  can  express  as  an  integer  linear  program.  Adding  additional  constraints  to 
the  master  problem  to  limit  the  political,  environmental  or  economic  impact  of  candidate 
delay  plans  [e.g.,  Reed,  1994]  is  straightforward,  as  is  allowing  more  than  one  type  of 
interdiction  action  for  each  task.  The  subproblem  could  be  modified  to  permit  a  single 
interdiction  action  to  delay  a  group  of  tasks  rather  than  only  one.  Our  subproblem  as¬ 
sumes  a  linear  cost  function  for  crashing,  although  it  could  be  modified  to  support  non¬ 
linear  cost  functions  [e.g.,  Vanhoucke,  et  al,  2002]  as  well. 
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There  is  no  requirement  that  the  master  problem  strietly  interdiet  the  proliferation 
subproblem.  We  ean  use  any  integer  linear  master  problem  we  please  to  manipulate  the 
proliferation  optimization  to  suit  our  needs.  For  instanee,  we  might  state  a  master  prob¬ 
lem  that  reschedules  a  proliferation  effort  to  maximize  its  eomplianee  with  telltale 
intelligenee  about  a  eovert  proliferator’s  aetions — evidenee  of  new  eonstruetion,  equip¬ 
ment  requisitions,  allianees  with  other  proliferators,  radiologieal  emissions,  power 
eonsumption  patterns,  eommunieations  intereepts,  et  cetera.  By  systematieally  enumerat¬ 
ing  alternate  proliferation  plans,  we  would  hope  to  identify  those  that  are  maximally 
eonsistent  with  what  we  observe.  This  immediately  yields  a  tool  for  evaluating  the  likely 
state  of  a  eovert  weapons  program  given  partial  information. 

Further  refinement  of  the  software  developed  here  eould  provide  valuable  fune- 
tionality.  Automating  the  ereation  of  visual  displays  sueh  as  Figures  1 1  and  13  would  be 
useful.  A  tool  to  visualize  ehanges  in  the  deeision  plan  in  response  to  ehanging  input  pa¬ 
rameters  would  be  helpful  as  well.  Migrating  the  optimization  model  from  GAMS  to  a 
general-purpose  programming  language  with  a  direetly  eallable  solver  paekage  sueh  as 
the  XA  eallable  subroutine  library  would  signifieantly  speed  solution  times  for  larger 
problems,  enabling  a  more  thorough  exploration  of  the  problem  spaee  for  the  deeision 
maker. 

We  are  suggesting  how  to  organize  and  manage  intelligenee  data  in  a  eohesive 
fashion  that  offers  visualization,  data  base  funetionality  and  interoperability  with  Miero- 
soft  Offiee.  The  optimization  is  new,  but  eleanly  hidden  from  view.  An  analyst,  or  even 
a  poliey  maker,  ean  pose  questions,  get  answers  and  gain  deep  insights  without  any 
mathematical  background  at  all. 
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LIST  OF  ABBREVIATIONS  AND  ACRONYMS 


AON 

Activity  on  node 

CPM 

Critical  path  method 

GAMS 

Generalized  Algebraic  Modeling  System 

HEU 

Highly  enriched  uranium 

IAEA 

International  Atomic  Energy  Agency 

MS  Project 

Microsoft  Project  software 

NERAC 

Nuclear  Energy  Research  Advisory  Committee 

NPT 

The  Treaty  on  the  Non-Proliferation  of  Nuclear  Weapons 

OTA 

Office  of  Technology  Assessment 

PERT 

Program  evaluation  and  review  technique 

SWU 

Separative  work  unit 

VBA 

Visual  Basic  for  Applications 

41 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


42 


LIST  OF  REFERENCES 


Ahuja,  R.  K.,  T.  L.  Magnanti  and  J.  B.  Orlin  [1993]  Network  Flows:  Theory,  Algorithms, 
and  Applications,  Prentice-Hall,  Inc. 

Brooke,  A.,  D.  Kendrick,  A.  Meeraus  and  R.  Raman  [1998]  GAMS:  A  User’s  Guide, 
GAMS  Development  Corporation,  Washington,  DC,  September. 

Crowston,  W.  and  G.  L.  Thompson  [1967]  “Decision  CPM;  A  Method  for  Simultaneous 
Planning,  Scheduling,  and  Control  of  Projects,”  Operations  Research,  Vol.  15,  No.  3, 
May- June,  pp.  407-426. 

Federation  of  American  Scientists  [2000]  “Osiraq  -  Iraq  Special  Weapons  Facilities,” 
[http;//www.fas.org/nuke/guide/iraq/facihty/osiraq.htm].  Accessed  January  2004. 

Gantt,  H.  L.  [1903]  “A  graphical  daily  balance  in  manufacture,”  ASME  Transactions, 
Vol.  24,  pp.  1322-1336. 

Hamming,  R.  W.  [1986]  Coding  and  Information  Theory,  2“^^  ed.,  Prentice-Hall,  p.  189. 

Harney,  R.  C.  [2003]  “Nuclear  Weapon  Development  Scenario,”  personal  communica¬ 
tion,  Naval  Postgraduate  School,  Monterey,  California,  December. 

International  Business  Machines  [1968]  “Project  Management  System,  Application  De¬ 
scription  Manual,”  H20-0210,  International  Business  Machines. 

Israeli,  E.  and  R.  K.  Wood  [2002]  “Shortest-Path  Network  Interdiction,”  Networks,  Vol. 
40(2),  pp.  97-111. 

Malcolm,  D.  G.,  J.  H.  Roseboom,  C.  E.  Clark  and  W.  Eazar  [1959]  “Application  of  a 
Technique  for  Research  and  Development  Program  Evaluation,”  Operations  Research, 
Vol.  7,  No.  5,  September-October,  pp.  646-669. 

Microsoft  Corporation  [2003a]  “Microsoft  Office,”  [http;//office.microsoft.com/home/ 
office. aspx?assetid=EX01085795],  Accessed  Eebruary  2004. 

Microsoft  Corporation  [2003b]  “Microsoft  Visual  Basic  for  Applications  Home  Page,” 
[http;//msdn.microsoft.com/vba/],  September  3,  Accessed  Eebruary  2004. 

Microsoft  Corporation  [2004]  “How  to  Determine  When  a  Shelled  Process  Ends,” 
[http://support.microsoft.com/default.aspx?scid=kb;en-us;178116&Product=vba],  Ac¬ 
cessed  January. 

Moder,  J.  J.,  C.  R.  Phillips  and  E.  W.  Davis  [1983]  Project  Management  with  CPM, 
PERT  and  Precedence  Diagramming,  3'^'*  ed..  Van  Nostrand  Reinhold  Company  Inc. 


43 


Nonproliferation  Policy  Education  Center  [2004]  “Iran;  Breaking  Out  Without  Quite 
Breaking  the  Rules?”  [http;//www.npec-web.org/projects/iranswu2.htm],  Accessed 
March  2004. 

Nuclear  Energy  Research  Advisory  Committee  (NERAC)  [2001]  “Technological  Oppor¬ 
tunities  To  Increase  the  Proliferation  Resistance  of  Global  Civilian  Nuclear  Power 
Systems  (TOPS),”  TOPS  Task  Eorce  Report,  January. 

Reed,  B.  K.  [1994]  Models  for  Proliferation  Interdiction  Response  Analysis,  MS  Thesis 
in  Operations  Research,  Naval  Postgraduate  School,  Monterey,  California,  September. 

Rhodes,  R.  [1986]  The  Making  of  the  Atomic  Bomb,  Simon  &  Schuster  Inc.,  p.  561. 

Rhodes,  R.  [1995]  Dark  Sun:  The  Making  of  the  Hydrogen  Bomb,  Simon  &  Schuster 
Inc.,  pp.  304-5. 


Spears,  D.  (ed.)  [2001]  “Technology  R&D  for  Arms  Control,”  Arms  Control  and  Non¬ 
proliferation  Technologies,  Eawrence  Livermore  National  Laboratory,  Spring,  p.  D. 

Sunset  Software  Technology  [2003]  “Sunset  Software  Technology  -  Optimization 
Tools,”  [http;//www.sunsetsoft.com/],  August  25,  Accessed  Eebruary  2004. 

United  Nations  [1968]  The  Treaty  on  the  Non-Proliferation  of  Nuclear  Weapons. 

United  Nations  [2004]  “Non-Proliferation  Treaty,”  [http;//disarmament.un.org;8080/ 
wmd/npt/].  Accessed  January. 

U.S.  Congress,  Office  of  Technology  Assessment  (OTA)  [1993]  “Technologies  Underly¬ 
ing  Weapons  of  Mass  Destruction,”  OTA-BP-ISC-115,  U.S.  Government  Printing  Office, 
December. 

Vanhoucke,  M.,  E.  Demeulemeester  and  W.  Herroele  [2002]  “Discrete  time/cost  trade¬ 
offs  in  project  scheduling  with  time-switch  constraints,”  Journal  of  the  Operational  Re¬ 
search  Society,  Vol.  53,  pp.  741-751. 

VBA-Programmer.com  [2004]  “Sorting  a  Listbox,”  [http;//www.vba-programmer.com/ 
VB_Code/Sortmg_Listbox.txt],  Accessed  January. 


44 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Teehnieal  Information  Center 

Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  Sehool 
Monterey,  California 

3.  Distinguished  Professor  Gerald  G.  Brown 

Naval  Postgraduate  Sehool 
Monterey,  California 

4.  Professor  R.  Kevin  Wood 
Naval  Postgraduate  Sehool 
Monterey,  California 

5.  Assoeiate  Professor  Robert  C.  Harney 
Naval  Postgraduate  Sehool 
Monterey,  California 

6.  Assoeiate  Professor  W.  Matthew  Carlyle 
Naval  Postgraduate  Sehool 
Monterey,  California 

7.  Anton  Rowe 

Naval  Postgraduate  Sehool 
Monterey,  California 

8.  LT  Erie  M.  Skroeh,  SC,  USN 
Marina,  California 


45 


